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The main goal of the paper is to give a short review on neutrino electromagnetic 
properties. In the introductory part of the paper a summary on what we really know 
about neutrinos is given: we discuss the basics of neutrino mass and mixing as well as the 
phenomenology of neutrino oscillations. This is important for the following discussion on 
neutrino electromagnetic properties that starts with a derivation of the neutrino electro- 
magnetic vertex function in the most general form, that follows from the requirement of 
Lorentz invariance, for both the Dirac and Majorana cases. Then, the problem of the 
neutrino form factors definition and calculation within gauge models is considered. In 
particular, we discuss the neutrino electric charge form factor and charge radius, dipole 
magnetic and electric and anapole form factors. Available experimental constraints on 
neutrino electromagnetic properties are also discussed, and the recently obtained experi- 
mental limits on neutrino magnetic moments are reviewed. The most important neutrino 
electromagnetic processes involving a direct neutrino coupling with photons (such as 
neutrino radiative decay, neutrino Cherenkov radiation, spin light of neutrino and plas- 
mon decay into neutrino-antineutrino pair in media) and neutrino resonant spin-flavor 
precession in a magnetic field are discussed at the end of the paper. 

PACS: 14.60. St, 13.15.+g 

1 Introduction 

The neutrino is a very fascinating particle which has remained under the focus of 
intensive investigations, both theoretical an experimental, for a couple of decades. These 
studies have given evidence of an ultimate relation between the knowledge of neutrino 
properties and the understanding of the fundamentals of particle physics. The birth of the 
neutrino was due to an attempt, by W. Pauli in 1930, to explain the continuous spectrum 
of beta-particles through "a way out for saving the law of conservation of energy" [T|. 
This new particle, called at first the "neutron" and then renamed the "neutrino", was 
an essential part of the first model of weak interactions (E. Fermi, 1934). Further im- 
portant milestones of particle physics, such as parity nonconservation (T.D. Lee, C.N. 
Yang and L. Landau, 1956) and the V — A model of local weak interactions (E. Sudar- 
shan, R. Marshak, 1956; R. Feynman, M. Gell-Mann, 1958), as well as the structure 
of the Glashow-Weinberg-Salam standard model, were based on the clarification of the 



* e-mail: giunti@to.infn.it, **e-mail: studenik@srd.sinp.msu.ru 



1 



specific properties of the neutrino. It has happened more than once that a novel dis- 
covery in neutrino physics stimulates far-reaching consequences in the theory of particle 
interactions. 

The neutrino plays a crucial role in particle physics because it is a "tiny" particle. 
Indeed, the scale of neutrino mass is much lower than that of the charged fermions 
(m Uf << irif, f = e, /i,r). The weak and electromagnetic interactions of neutrinos 
with other particles are really very weak. That is a reason for the neutrino to fall under 
the focus of researchers during the latest stages of a particular particle physics evolu- 
tion paradigm when all of the "principal" phenomena have been already observed and 
theoretically described. 

Neutrino electromagnetic properties, that is the main subject of this paper, are of 
particular importance because they provide a kind of bridge to "new physics" beyond 
the standard model. In spite of reasonable efforts in studies of neutrino electromagnetic 
properties, up to now there is no experimental confirmation in favour of nonvanishing 
neutrino electromagnetic characteristics. The available experimental data in the field 
do not rule out the possibility that neutrinos have "zero" electromagnetic properties. 
However, in the course of the recent development of knowledge on neutrino mixing and 
oscillations, supported by the discovery of flavor conversions of neutrinos from different 
sources, non-trivial neutrino electromagnetic properties seem to be very plausible. 

The structure of the paper is as follows. In the first part of the paper we summa- 
rize what we really know about neutrinos: the basics of neutrinos mass and mixing are 
discussed, as well as the phenomenology of neutrino oscillations. This introductory part 
is important for understanding the second part of the paper that is devoted to electro- 
magnetic (in the sense mentioned above still "unknown") properties of neutrinos. We 
start discussing neutrino electromagnetic properties deriving the neutrino electromag- 
netic vertex function in the most general form for both the Dirac and Majorana cases. 
Then, we consider the neutrino electric charge form factor and charge radius, magnetic, 
electric and anapole form factors. We discuss the relevant theoretical items as well as 
available experimental constraints. In particular, the neutrino magnetic and electric mo- 
ments, in both theoretical and experimental aspects, are discussed in detail. In Section [31 
the most important neutrino electromagnetic processes involving the direct neutrino cou- 
plings with photons (such as neutrino radiative decay, neutrino Cherenkov radiation, 
spin light of neutrino and plasmon decay into neutrino-antineutrino pair in media) and 
neutrino resonant spin-flavor precession in a magnetic field are discussed. 

2 What we know about neutrino 

In the Standard Model of electroweak interaction, forged in the 60's by Glashow, 
Weinberg and Salam [2H3], neutrinos are massless by construction (see Ref. [5]). This 
requirement was motivated by the low experimental upper limit on the neutrino mass 
(see Fig. [TJ and by the theoretical description of neutrinos through massless left-handed 
Weyl spinors in the two-component theory of Landau, Lee and Yang, and Salam [SHE], 
which prompted the V — A theory of charged- current weak interactions of Feynman and 
Gell-Mann, Sudarshan and Marshak, and Sakurai [9HTTj . 

The massless of neutrinos in the Standard Model is due to the absence of right-handed 
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Figure 1: Order of magnitude of the masses of leptons and quarks. 



Table 1: Eigenvalues of the weak isospin /, of its third component ^3, of the hypercharge 
Y, and of the charge Q = I3 + Y/2 of the lepton and Higgs doublets and singlets in the 
extension of the Standard Model with the introduction of right-handed neutrinos. 
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neutrinos, without which it is not possible to have Dirac mass terms, and to the absence 
of Higgs triplets, without which it is not possible to have Majorana mass terms. In the 
following we will consider the extension of the Standard Model with the introduction 
of three right-handed neutrinos. We will see that this seemingly innocent addition has 
the very powerful effect of allowing not only Dirac mass terms, but also Majorana mass 
terms for the right-handed neutrinos which induce Majorana masses for the observable 
neutrinos. 

Table [1] shows the values of the weak isospin, hypercharge, and electric charge of the 
lepton and Higgs doublets and singlets in the extended Standard Model under consider- 
ation. For simplicity, we work in the flavor base in which the mass matrix of the charged 
leptons is diagonal. Hence, e, /1, r are the physical charged leptons with definite masses. 
The three singlet neutrinos are often called sterile^ since they do not take part in weak 
interactions, in contrast with the standard active neutrinos v e , v^, v T . 

2.1 Dirac mass term 

The fields in Tab. [1] allow us to construct the Yukawa Lagrangian term 

3 

£y = - $^fc^$^ fc R + H.c., (1) 

a=e,fj,,r k=l 
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where Y is a matrix of Yukawa couplings and $ = zo"2$*. In the Standard Model, a 
nonzero vacuum expectation value of the Higgs doublet, 

1 (6 



induces the spontaneous symmetry breaking of the Standard Model symmetries SU(2) L x 
U(l) y — > U(1)q. In the unitary gauge, the Higgs doublet is given by 

= "4 („ , J) > (3) 



^2 + 

where H(x) is the physical Higgs field. From the Yukawa Lagrangian term in Eq. (CD) we 
obtain the neutrino Dirac mass term 

3 

-2b = — y= J2 YakW ^ + H - c - • ( 4 ) 

* a=e,[i,T k=l 

Since the matrix Y is, in general, a complex 3x3 matrix, the flavor neutrino fields u e , 
u^, v T do not have a definite mass. The massive neutrino fields are obtained through the 
diagonalization of j£fo. This is achieved through the transformations 



3 

V a L = U a k VkL , V Sj r = Vjk VkR , (5) 

k=l k=l 



with unitary matrices U and V which perform the biunitary diagonalization 



-j=(rfYV) kj =m k 6 hj , (6) 
with real and positive masses m k . The resulting diagonal Dirac mass term is 

3 3 
=^b = — m k ~VkL v kR + H.c. = — VTL k T% Vj. , (7) 

fc=l fe=l 
with the Dirac fields of massive neutrinos v k — v kL + VkR- 



2.2 Dirac— Majorana mass term 

In the above derivation of Dirac neutrino masses we have implicitly assumed that the 
total lepton number is conserved. If this assumption is lifted, neutrino masses receive 
an important contribution from the Majorana mass term of the right-handed singlet 
neutrinos, 

1 3 

k,j=i 

where C is the charge-conjugation matrix. The mass matrix M R is complex and symmet- 
ric. 
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The Majorana mass term in Eq. (jHJ) is allowed by the symmetries of the Standard 
Model, since right-handed neutrino fields are invariant. On the other hand, an analogous 
Majorana mass term of the left-handed neutrinos, 

*l=\ E ^M^ + H.c., (9) 

o,/3=e,//,T 

is forbidden, since it has I3 = 1 and Y = —2. There is no Higgs triplet in the Standard 
Model to compensate these quantum numbers. 

In the extension of the Standard Model with the introduction of right-handed neutri- 
nos, the neutrino masses and mixing are given by the Dirac-Majorana mass term 

=2d+m = + ■ (10) 

The neutrino fields with definite masses are obtained through the diagonalization of 
J^d+m- It is convenient to define the vector Nl of 6 left-handed fields 

Nl = (u eL , v^l, u tL , (u siR ) c , (u S2 r) c , (v S3 r) c ) , (11) 

with the charge-conjugated sterile neutrino fields {v s r) c = CVJr 1 . The Dirac-Majorana 
mass term in Eq. ffTO]) can be written in the compact form 

^d+m = \n T l & M d+m N l + H.c. , (12) 
with the N x N symmetric mass matrix 

M ° +M s (A M« ) • < 13 > 

where 

M D = -^r. (14) 



Notice that the Dirac-Majorana mass term in Eq. (jl2j) has the structure of a Majorana 
mass term. Therefore, it will not be a surprise to find in the following that the 6 massive 
neutrinos obtained from the diagonalization of J2?d+m are Majorana particles. 
Equation ( TT2"j) is diagonalized through the unitary transformation 

N L = Vn L , with nl = (u 1L , u 6L ) . (15) 

The unitary matrix V is chosen in order to diagonalize the symmetric mass matrix M D+M : 

yT M D+M y = M ^ where M ^ = ^ (k,j=l,...,6), (16) 

with real and positive masses m&. In this way, the Dirac-Majorana mass term in Eq. ffT2|) 
can be written in terms of the massive fields as 

1 1 6 1 6 

-2d+m = -n T L C ] Mn L + H.c. = - ^ m k v T kL C ] u kL + H.c. = - ^ m k v T k C ] v k , (17) 

fc=l k=l 

where v k = v k L + v£ L are Majorana fields which satisfy the constraint u% = v k . Hence, 
a general result of the diagonalization of a Dirac-Majorana mass term is that massive 
neutrinos are Majorana particles. 

Note that in the limit M R = we recover the Dirac case, since there are three pairs of 
degenerate mass eigenvalues. Each pair of massive Majorana fields with the same mass 
corresponds to a Dirac field (see Ref. [5]). For example, if the massive Majorana fields v\ 
and z/ 4 have the same mass mi, they correspond to the Dirac field (v\ + iv^)/y/2. 
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2.3 Weak interactions 

The mixing of neutrinos in Eq. (jl5[) is observable through its effect in weak interac- 
tions. Let us first consider the leptonic weak charged current 



Hence, only the rectangular submatrix of V composed by the first three rows is relevant for 
charged-current weak interactions. Each of the 6 massive neutrinos partake in charged- 
current weak interactions if the elements the first three rows and corresponding column 
of the mixing matrix are not negligibly small. In Section 12.41 we will see that in the 
celebrated see-saw mechanism the effective number of light Majorana massive neutrinos 
which take part in weak interactions is reduced to three. 
The neutrino weak neutral current is given by 



Unless the mixing of 6 neutrinos reduces to an effective mixing of 3 neutrinos, as in 
the see-saw mechanism discussed in Section I2.4[ J2 a=e „ T V* k V a j ^ 5kj and there can be 
neutral-current transitions among different massive neutrinos (no GIM mechanism [12])- 

2.4 See-saw mechanism 

The order of magnitude of the elements of the Dirac mass matrix M D in Eq. (fT4|) 
is expected to be smaller than v ~ 10 2 GeV, since the Yukawa couplings are expected 
not to be unnaturally large. In general, since a Dirac mass term is forbidden by the 
symmetries of the Standard Model, it can arise only as a consequence of symmetry 
breaking and Dirac masses are proportional to the symmetry-breaking scale. This fact 
is often summarized by saying that Dirac masses are protected by the symmetries of 
the Standard Model. On the other hand, since the Majorana mass term in Eq. (jSJ) is a 
Standard Model singlet, the elements of the Majorana mass matrix Mr are not protected 
by the Standard Model symmetries. It is plausible that the Majorana mass term Jz?# 
is generated by new physics beyond the Standard Model and the right-handed chiral 
neutrino fields v s r belong to nontrivial multiplets of the symmetries of the high-energy 
theory. In this case, the elements of the mass matrix Mr are protected by the symmetries 
of the high-energy theory and their order of magnitude corresponds to the breaking scale 
of these symmetries, which may be as large as the grand unification scale, of the order of 
10 14 -10 16 GeV. The mass matrix can be diagonalized by blocks, up to corrections of the 
order (M^M : 
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The light 3x3 mass matrix Mu g ht and the heavy 3x3 mass matrix Mh ea vy are given by 
Might — —M dT (M R )~ l M D , M hcavy ~M R . (22) 

The heavy masses are given by the eigenvalues of M R , whereas the light masses are given 
by the eigenvalues of M light , whose elements are suppressed with respect to the elements 
of the Dirac mass matrix M D by the very small matrix factor M D (M R )~ l . This is the 
celebrated see-saw mechanism, which explains naturally the smallness of light neutrino 
masses (see Fig. [T]). Notice, however, that the values of the light neutrino masses and 
their relative sizes can vary over wide ranges, depending on the specific values of the 
elements of M D and M R . 

Since the off-diagonal block elements of W are very small, the three flavor neutrinos 
are mainly composed by the three light neutrinos. Therefore, the see-saw mechanism im- 
plies the effective low-energy mixing of three Majorana neutrinos with an approximately 
unitary 3x3 mixing matrix U composed by the first 3 rows and the first 3 columns of V. 

2.5 Three-neutrino mixing 

In the case of three-neutrino mixing, the three left-handed flavor neutrino fields u e £, 
VfiL-, v tL which partake in weak interactions are unitary linear combinations of three 
left-handed massive neutrino fields vu, ^2L, ^3L- 

3 

VaL = ^ Uak VkL = e ' ^' T ) ' ( 23 ) 

k=l 

where U is a 3 x 3 unitary mixing matrix. Motivated by the see-saw mechanism, we 
consider Majorana massive neutrinos. The deviation from the unitarity of U in the see- 
saw mechanism is negligible. In this approximation, the expression in Eq. (ITS]) for the 
leptonic weak charged current reduces to 

3 

f cc = 2 ^TlY* l = 2 ^VkZUrt-faL, (24) 

a=e,n,T a=e,fi,r k=l 

and the expression in Eq. (Tl9|) for the leptonic weak neutral current reduces to 

3 

■7nc = ^2 v *l1 Pv <*l= ^VkEYvjL, (25) 

a=e,fi,T k,j=l 

where the unitarity of U implies the absence of neutral-current transitions among different 
massive neutrinos (GIM mechanism |12j). 

Notice that in the case of three-neutrino mixing the mixing matrix U enters only 
in the leptonic weak charged current j^ c . Hence, it is observable only through weak 
charged-current interactions. 

The unitary matrix U can be parameterized in terms of 3 mixing angles and 6 phases. 
However, 3 phases are unphysical, because they can be eliminated by rephasing the three 
charged lepton fields in j^ c . In the case of Majorana massive neutrinos, no additional 
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phase can be eliminated, because a Majorana mass term vj.& Vk is not invariant under 
rephasing of v^. On the other hand, in the case of Dirac massive neutrinos, two additional 
phases can be eliminated by rephasing the massive neutrino fields. Hence, the mixing 
matrix has 3 physical phases in the case of Majorana massive neutrinos or 1 physical 
phase in the case of Dirac massive neutrinos. In general, in the case of Majorana massive 
neutrinos U can be written as 

U = U D D M , (26) 

where U D is a Dirac unitary mixing matrix which can be parameterized in terms of three 
mixing angles and one physical phase, called Dirac phase, and D M is a diagonal unitary 
matrix with two physical phases, usually called Majorana phases. 
The standard parameterization of U D is 

/ ci 2 c 13 S12C13 si 3 e- 4,5l3 \ 

U D = -s 12 c 23 - c 12 s 23 s 13 e lSl3 c 12 c 23 - s 12 s 23 s 13 e lSl3 s 23 c 13 , (27) 
V s 12 s 23 - c 12 c 23 s 13 e tSl3 -c 12 s 23 - si 2 c 23 s 13 e lSl3 c 23 c 13 J 

where c a b = cos$ a j, and s a b = sm$ ab . $ 12 , $13, $23 are the three mixing angles (0 < $ a b < 
7r/2) and Si 3 is the Dirac phase (0 < 5\ 3 < 2ir). 

The diagonal unitary matrix D M can be written as 

D M = diag(e iAl , e iX2 , e iXs ) , with A x = . (28) 

The phases X 2 and A3 are the two physical Majorana CP-violating phases. Since all 
measurable quantities depend only on the differences of the three phases Ai, A2, A3, the 
choice Ai = is a matter of convention and other choices are equivalent from the physical 
point of view. In fact, rephasing all the charged lepton fields in by e llf , we have 
e iAfe — > e t( - Xk ~'~ p \ whereas e*^ - ^ remains constant. 

All the phases in the mixing matrix violate the CP symmetry. Since the Majorana 
phases are observable only in processes which are allowed only in the case of Majorana 
neutrinos, as neutrinoless double-/? decay, in most observable processes CP violation is 
generated by the Dirac phase. The size of this CP violation can be quantified in a 
parameterization-invariant way by the Jarlskog invariant [T3"HT%] 

J = lm(U e2 U, 3 C/ e * 3 U; 2 ) = Im(f/ e D 2 U% U% U%) . (29) 

The second equality is due to the invariance of J under the rephasing U a k — > e~ %lfa U a ke l ^ k , 
with arbitrary phases <p a and ipk, which implies that the Majorana phases do not con- 
tribute. Using the unitarity of the mixing matrix, it can be shown that 

J am = lm(U ak U Pj U*j U; k ) =±J. (30) 

Therefore, in the case of three-neutrino mixing | J\ quantifies Dirac CP violation inde- 
pendently from the parameterization of the mixing matrix. In the parameterization in 
Eq. (I27p . we have 

J = Ci2Si2C2 3 S23C?3Si 3 sin5i3 = - sin 2^12 sin 2^ 2 3 cos $13 sin 2$i3 sin5i 3 . (31) 
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2.6 Neutrino oscillations 



Neutrino oscillations was proposed by B. Pontecorvo in the late 1950s in analogy 
with K°-K° oscillations [TjJl ED] The oscillations are generated by the interference of 
the phases of different massive neutrinos, which are produced and detected coherently 
because of their very small mass differences. 

Let us consider a flavor neutrino state 



Wa) = ^Kk k) 



(32) 



which describes a neutrino with flavor a created in a charged- current weak interaction 
process from a charged lepton a~ or together with a charged antilepton a + (a = e, fi, r). 
The presence of the weight U* k for \v k ) in the flavor state \u a ) is due to the decomposition 
in Eq. (j23J) of the leptonic charged current Jq C in terms of the massive neutrino contribu- 
tions, which contain the creation operators of massive neutrinos. Additional coefficients 
due to different effects of neutrino masses in the interaction process are negligible in 
neutrino oscillation experiments. 

The massive neutrino states \v k ) are eigenstates of the free Hamiltonian with energy 
eigenvalues 

E k = \/\Pk\ 2 + ml, (33) 



where p k is the respective momentum. Since the massive neutrinos evolve in space-time 
as plane waves, the space-time evolution of the flavor state in Eq. (1321) is 



\v a {L 1 T)) = Y J K k e~ lEkT+iPk - L H)= J2 [T, U ^e- lEkT+iPkL U, k )\u,), (34) 

k P=e,fi,T \ k / 

where we have used the unitarity of the mixing matrix for inverting the relation in 
Eq. ( 13"2"j) . One can see that the phase differences of different massive neutrinos gener- 
ate flavor transitions with probability 



P Va ^„JL,T) = \(v f) \v a (L,T))f 



(35) 



Since the source-detector distance L = |X| is macroscopic, we can consider all massive 
neutrino momenta p k aligned along L. Moreover, taking into account the smallness of 
neutrino masses, in oscillation experiments in which the neutrino propagation time T 
is not measured it is possible to approximate T = L. With these approximations, the 
phases in Eq. fl35|) reduce to 



E k T + p k L = - (E k - p k ) L 



El 



Pk 



E k +p k 



L 



mi 



E k +p k 



L 



(36) 



at lowest order in the neutrino masses. Here, p k = |pk| and E is the neutrino energy 
neglecting mass contributions. Equation ( 136|) shows that the phases of massive neutrinos 
relevant for oscillations are independent of the values of the energies and momenta of 
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different massive neutrinos, because of the relativistic dispersion relation in Eq. f l33|) . 
The flavor transition probabilities are 



P„ a ^{L, E) = 8 a p-±J2 MKk Upu U aj U* pj ) sin 2 ( 



Hi 

AE 



k>j 

x ^ /Am 2 L\ 

-2^2j afikj sml-^-\ . (37) 

k>j v ' 

The C and T conjugated flavor transition probabilities are given by 

The survival probabilities (a = (3) are CP- invariant (a consequence of CPT symmetry), 

Pv a ^V a = Ppa^tUa ) (39) 

whereas CP violation is observable in flavor transitions by measuring the asymmetries 

4* = P» a -+u, - Pp.-*, («^). (40) 

CPT symmetry implies that the CP asymmetries are equal to the corresponding T asym- 
metries: A£p = —A^p = A^, with A^p = P Va ^ Vf) — P Vfj ^ Va and A^ = Pp a ^p p — Pp p -+p a - 
In the approximation of two-neutrino mixing, in which one of the three massive neu- 
trino components of two flavor neutrinos is neglected, the mixing matrix reduces to 

TT ( cos i? sin$\ , . 

u = n o , 41 

sin 17 cos v J 

where $ is the mixing angle (0 < d < n/2). In this approximation, there is only one 
squared-mass difference Am 2 and the transition probability is given by 



Am 2 L 
AE 



P Va -> Vp (L, E) = sin 2 2§ sin 2 ( ) (o^). (42) 



In the case a = (3, the survival probability is 



P» a ^ a (L, E) = l- sin 2 2tf sin 2 (^pj • (43) 

These simple expressions are often used in the analysis of experimental data. 

When neutrinos propagate in matter, the potential generated by the coherent forward 
elastic scattering with the particles in the medium (electrons and nucleons) modifies 
mixing and oscillations [21]. In a medium with varying density it is possible to have 
resonant flavor transitions [22]. This is the famous MSW effect. 

The effective potentials for u a and v a are, respectively, 

V a = V C c5 ae + V NC , V a =-V a , (44) 
with the charged-current and neutral-current potentials 



V2G F N ej v NC = -\V2G F N nj (45) 
10 
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Figure 2: Feynman diagrams of the coherent forward elastic scattering processes that 
generate the CC potential Vcc through W exchange and the NC potential Vnc through 
Z exchange. 

generated by the Feynman diagrams in Fig. [2J Here N e and N n are the electron and 
neutron number densities in the medium (in an electrically neutral medium the neutral- 
current potentials of protons and electrons cancel each other). In normal matter, these 
potentials are very small, because 

v^G^7.63xlO- 14 ^^, (46) 

where N\ is Avogadro's number. 

Let us consider, for simplicity, two-neutrino v^-v^ mixing. In general, a neutrino 
produced at x = is described at a distance a; by a state 

\v(x)) = (p e {x) \v e ) + if^x) |r/ M ) . (47) 

The evolution of the flavor amplitudes <p e (x) and (p^{x) with the distance x is given by 
the differential equation [21] 



d_ (<p e (x)\ = f^sm 2 $ + V e ^sin2tf \ ( Ve {x) . 
dx\<PM) V ^sin2tf ^ C0S 2# + vJ V^»/ ' 



For an initial u e , the boundary condition for the solution of the differential equation is 

(49) 



W°)/ v° 



and the probabilities of v e — \ transitions and v e survival are, respectively, 

P^(x) = |^(x)| 2 , P v ^ VB (x) = \<p e (x)\ 2 = 1 - P^(x) . (50) 

The evolution equation (f48l) has the structure of a Schrodinger equation with the 
effective Hamiltonian matrix 

Am 2 1 1 /- Am 2 cos 2tf + 2EV C c Am 2 sin 2$ \ , . 

H ~^E~ + 2 CC+ NC + 4E\ Am 2 sin 2tf Am 2 cos 2d - 2EV CC ) ' ( j 
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This matrix can be diagonalized by the orthogonal transformation 



Am 2 1 1 
U^UU M = — + - V cc + Vkc + |£ diag(- Am^, Am 2 M ) . (52) 



The orthogonal matrix 

sin i? M cos t?M 
is the effective mixing matrix in matter, and 



/ cos^m sin^M \ ,k Q \ 



Am^ = y (Am 2 cos 2d - 2i?Vcc) 2 + (Am 2 sin 2d) 2 (54) 

is the effective squared-mass difference. The effective mixing angle in matter is given 
by 

tan 2$ 

tan2# M = ttwtf • (55) 



2EV ( 



cc 



Am 2 cos 2$ 

The most interesting characteristic of this expression is that there is a resonance [22] 
when 

Am 2 

^cc = ^cos2tf, (56) 

which corresponds to the electron number density 

^ R Am 2 cos2t? 
e 2V2EG F 

At the resonance the effective mixing angle is equal to 7r/4, i.e. the mixing is maximal, 
leading to the possibility of total transitions between the two flavors if the resonance 
region is wide enough. 

In general, the evolution equation ( l4"g|) must be solved numerically or with appropriate 
approximations. In a constant matter density, it is easy to derive an analytic solution, 
leading to the transition probability 

iW*) = sin 2 2d M sin 2 (^§^) , (58) 

which has the same structure as the two-neutrino transition probability in vacuum in 
Eq. (1421) . with the mixing angle and the squared-mass difference replaced by their effective 
values in matter. 



2.7 Phenomenology 

Since the 1960s it has been known, mainly through the insight of Pontecorvo [23] , that 
neutrino oscillations can be revealed not only in terrestrial neutrino experiments, but also 
in experiments which are sensitive to neutrinos coming from astrophysical sources. The 
largest astrophysical neutrino flux on Earth coming from the Sun has been measured by 
several experiments, starting with the pioneering Homestake experiment [21] , which first 
observed the deficit of electron solar neutrinos with respect to the Standard Solar Model 
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Figure 3: [28] Ratio of the background and geoneutrino-subtracted v e spectrum to the 
expectation for no-oscillation as a function of L /E. L is the effective baseline taken as 
a flux- weighted average (L = 180 km). 



prediction (see Ref. In 2002 the SNO experiment [2B] has shown that the solar 

neutrino problem is due to v e — > v T transitions. The reactor long-baseline KamLAND 
experiment established at the end of 2002 [27] that these v e — > v T transitions are 
due to neutrino oscillations. The evidence for oscillations obtained by the KamLAND 
experiment is illustrated in Fig. |3j The current solar and KamLAND data are fitted 
well by effective two- neutrino oscillations, including MSW [2"Tjl2"2"] effects of neutrino 
propagation in matter, with the solar squared-mass difference and mixing angle [28~| 

A m 2 UN = (7.59 ± 0.21) x 10~ 5 eV 2 , tan 2 $ S un = 0.47±g;gg . (59) 

The allowed regions in the tan 2 ^suN _ Amg UN plane obtained from KamLAND data and 
the data of solar neutrino experiments (Homestake [21], GALLEX/GNO [29], SAGE [30] , 
Super-Kamiokande [21], SNO [22], Borexino [22]) are shown in Fig. |H Figure E] shows the 
results of the analysis of solar neutrino data alone and in combination with KamLAND 
data. 

The atmospheric neutrino anomaly was discovered in the late 1980s in the Kamiokande 
[21] and 1MB [25] experiments. In 1998 the Super-Kamiokande experiment found a model 
independent evidence of muon (anti) neutrino disappearance in atmospheric neutrino data 
[5B] . Atmospheric neutrinos are produced by the decay of pions and muons created by 
the interactions of cosmic rays with the nuclei in the atmosphere. Since at energies 
higher than about 1 GeV the flux of atmospheric neutrinos is approximately isotropic, 
the corresponding number of events generated in a detector by atmospheric neutrinos 
must be the same in any direction. The Super-Kamiokande collaboration measured the 
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Figure 4: [28] Allowed regions for v e — )• z/ M , i/ T oscillation parameters from KamLAND 
and solar neutrino data. The side-panels show the Ax 2 -profiles for KamLAND (dashed) 
and solar experiments (dotted) individually, as well as the combination of the two (solid). 



up-down asymmetry of high-energy ^-induced events [36] 

A up-down = (^ZR) = -0.296 ± 0.049 , (60) 

where U and D are, respectively, the neutrino fluxes integrated in the ranges —1 < 
cos 9 Z < —0.2 and 0.2 < cos# z < 1 (8 Z is the angle between the zenith and the neutrino 
arrival direction). Since the measured asymmetry deviates from zero by about 6a, the 
model-independent evidence of an atmospheric neutrino anomaly is indisputable. The 
negative value of ^4^p- down indicates that muon (anti)neutrinos coming from the opposite 

hemisphere disappear, most likely because of ( z/^ — > v T oscillations with the mixing pa- 
rameters in Fig. |6j since atmospheric electron (anti)neutrinos do not show any anomalous 

behavior. This interpretation has been confirmed by the independent observations of 
disappearance in the accelerator long-baseline experiments K2K [38] and MINOS [39] 
which are generated by the same values of the mixing parameters, as shown in Fig. [71 
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Figure 5: [32] Allowed regions for v e — > v T 
lar+KamLAND (c) data. The best-fit points 
from solar data (b) and Am 2 = 7.59 x 10~ 5 
data (c). 
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Figure 6: [37] Allowed regions for'i^ — y y v T oscillation parameters from Super-Kamiokande 
atmospheric neutrino data obtained with a zenith angle analysis (left) and a L/E analysis 
(right). The best-fit values of the oscillation parameters are: sin 2 2$ = 1.02, Am 2 = 
2.1 x 10~ 3 eV 2 (zenith) and sin 2 2$ = 1.04, Am 2 = 2.2 x 10~ 3 eV 2 (L/E). 



(-) 



From MINOS data [39], 

Am 2 ATM = (2.43 ± 0.13) x 10" 3 eV 2 , sin 2 2^ A tm > 0.90 (90% C.L.) . (61) 

The different values of solar and atmospheric squared-mass differences in Eqs. (1591) and 
(ISTf imply that two-neutrino mixing, with one squared-mass difference, is not sufficient 
for the description of all oscillation data. Moreover, all three neutrino flavors are involved 
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Figure 7: [39] Allowed regions for ^-disappearance oscillation parameters from K2K 
and MINOS data. 
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Figure 8: The two three- neutrino schemes allowed by the hierarchy Amg OL <C Am| TM . 



in the observed oscillations (solar and reactor v e — > v^,v T ] atmospheric and accelerator 

v v ~~ v -n- Therefore we must consider the mixing of three neutrinos in Eq. fT2B"|) . The 
observed hierarchy Am| L <S Am\ TM can be accommodated in the normal and inverted 
three-neutrino mixing schemes shown schematically in Fig. |HJ We choose the arbitrary 
labeling numbers of the massive neutrinos in order to have Amg OL = Am^ and Am\ TM = 
lAmjjiJ, with Am^ <C Afflgj ~ Am| 2 . 

In principle, the analysis of neutrino oscillation data in a three-neutrino mixing 
framework could yield results which are different from those obtained in a two-neutrino 
mixing approximation. However, in practice the two-neutrino mixing approximation is 
quite accurate, because the only element of the mixing matrix which affects both so- 
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lar and atmospheric neutrino oscillations, U e3 , is small. This information comes from 
the results of the CHOOZ [10] and Palo Verde [4TJ reactor long-baseline experiments, 
which excluded v e disappearance generated by Am\ TM with the effective mixing angle 
sin 2 2^f = 4 \U e3 \ 2 (1 - \U e3 \ 2 ) = sin 2 2& 13 g2]. The results of a global analysis of neu- 
trino oscillation data are shown in Fig. |9j The 90% C.L. (3<r) bounds for sin 2 $ 13 are jl3] 

( 0.060 (0.089) (solar+KamLAND) , 
sin 2 ^i 3 << 0.027 (0.058) (CHOOZ+atm+K2K+MINOS) , (62) 
{ 0.035 (0.056) (global data) . 

Therefore, in practice we have 

$SOL ^ $12 , $ATM ^ $23 , (63) 

and the results in Eqs. (l59]) and flBTj) apply to three- neutrino mixing. 

So far we have considered only neutrino oscillation data, which give information on 
neutrino mixing and the differences of neutrino squared masses. The absolute scale of 
neutrino masses must be determined with other means. However, since we know the 
squared-mass differences from Eqs. (159]) and ( l6Tj) . it is possible to express the neutrino 
masses as functions of only one unknown parameter representing the absolute mass scale. 
Figure [TU] shows the values of the three neutrino masses as functions of the lightest mass, 
which is rrii m the normal scheme and m 3 in the inverted scheme. One can see that 
in both schemes there is quasidegeneracy of the three masses when nil ~ 777-2 — ^3 3> 
\J Am 2 YFM ~ 5 x 10~ 2 eV. In this case, it is very difficult to distinguish the two schemes. 
On the other hand, the two schemes have very different features if the lightest mass is 
much smaller than Am| TM . In this case, in the normal scheme there is a hierarchy of 
masses: m\ m 2 m 3 . In the inverted scheme there is a so-called inverted hierarchy 
7773 rrii — m-2 in which nil and 7712 are quasidegenerate. In fact, in the inverted scheme 
nil and 7772 are always quasidegenerate, because their separation is due to the small solar 
squared-mass difference A?77g OL . Let us note that, independently of the mass scheme, at 
least two neutrinos are massive, with masses larger than about 8 x 10 -3 eV. 
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Lightest Mass: mi [eV] Lightest Mass: 7713 [eV] 



Figure 10: Values of neutrino masses as functions of the lightest mass mi in the normal 
scheme and in the inverted scheme. Solid lines correspond to the best-fit values of 
Amg UN and Am\ TM . Dashed lines enclose 3a ranges. 

The most reliable method for the determination of the absolute value of neutrino 
masses is the kinematic measurement of neutrino masses in interactions. Currently, the 
best limit is obtained in tritium /3-decay experiments, which are sensitive to the effective 
mass 

= ^\U ek ?ml. (64) 

The current bound on mp was obtained in the Mainz [H] and Troitzk [15] experiments: 

m /3 <2.3eV (95%C.L.). (65) 

Figure fTTl shows the comparison of this bound with the possible value of mp in the normal 
and inverted schemes as a function of the lightest mass. 

Another very important process which is sensitive to the absolute scale of neutrino 
masses is neutrinoless double-/3-decay, which occurs only if massive neutrinos are Majo- 
rana particles. Neutrinoless double-/3-decay depends on the effective Majorana mass 

3 

m 2 p = U ek m k ■ (66) 
fc=i 

The best limit on m 2( g, obtained in the Heidelberg-Moscow 76 Ge experiment, [16] is 

Kal <0.3-1.0eV, (67) 

where the large uncertainty is of theoretical nuclear physics origin. Figure [12] shows 
the comparison of this bound with the possible value of m 2 /3 in the normal and inverted 
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Figure 11: Effective neutrino mass mp in tritium /3-decay experiments as a function of 
the lightest mass (mi in the normal scheme and in the inverted scheme; see Fig. [8]). 
Middle solid lines correspond to the best-fit values of Awig UN and Am\ TU . Extreme solid 
lines enclose 3a ranges. Dashed lines show the best-fit values and 3a ranges of individual 
masses. In the inverted scheme, the best-fit values and 3a ranges of mi and m 2 are 
practically the same and coincide with the best-fit value and 3a range of mp. 




Figure 12: Absolute value \mip\ of the effective Majorana neutrino mass in 2/?o I /-decay 
as a function of the lightest mass mi in the normal scheme (a) and in the inverted 
scheme (b). The two horizontal dotted lines correspond to the extremes of the upper 
bound range in Eq. ( 16Tj) . The two vertical dotted lines show the corresponding upper 
bounds for mi (a) and m 3 (b). 
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Figure 13: Neutrino electromagnetic vertex function. 

scheme as a function of the lightest mass. The unshaded strip within the shadowed bands 
can be obtained only in the case of CP violation. 

Let us finally mention that the evolution of the Universe depends on the values of 
neutrino masses (see Ref. [5]). Current cosmo logical data limit the sum of neutrino masses 

by m 

3 

J^m fc <0.2-0.7eV, (68) 
fc=i 

in the framework of the very successful flat ACDM model. 



3 Neutrino electromagnetic properties 

The importance of neutrino electromagnetic properties was first mentioned by Pauli 
just in 1930 when he postulated the existence of this particle and discussed the possibility 
that the neutrino might have a magnetic moment. Systematic theoretical studies of 
neutrino electromagnetic properties have started after it was shown that in the extended 
Standard Model with right-handed neutrinos the magnetic moment of a massive neutrino 
is, in general, nonvanishing and that its value is determined by the neutrino mass jHH53] . 

Neutrino electromagnetic properties are of particular importance because they are di- 
rectly connected to fundamentals of particle physics. For example, neutrino electromag- 
netic properties can be used to distinguish Dirac and Majorana neutrinos (see [5T |Iq^H5T] 
for the correspondent discussion) and also as a probe of new physics that might exist 
beyond the Standard Model (see, for instance, [58|l59]). 

Consider the matrix element of the electromagnetic current between the fermion initial 
state i>{p) and final state tp(p') can be presented in the form 

< W)\J* M mP) >= u(p')A„(qJ)u(p), (69) 

where = p'^ — p^, = p'^ + p M . The matrix element between the spinors of the elec- 
tromagnetic vertex function A^(q, I) (Fig. [TBI) should be a Lorentz vector (requirement of 
Lorentz-covariance). In constructing the covariant operator A A1 (g, /) we recall □ that there 
are 16 linearly independent traceless (with the exception of the unit matrix) matrices, 

1, 75, 7/,, 757m, <V, ( 70 ) 

Cfiu = f b7i,7v]- There are in addition also the metric tensor g^, two vectors q^ and 1^, 
and the anti-symmetric tensor e^ ual that can be used. 

1 A rather pedagogical discussion on the electromagnetic form factors of spin-^ particles is given in |60j . 
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There are three sets of operators from which A At (g, I) can be formed. In the first set 
the Lorentz index is carried by the vectors g M and / M , 

lip, 75^, 7 5 Z M . (71) 

There is another set of the same type, 

75?/., Ib/iq^ 7s v a pq a l^q^ (72) 

and the correspondent terms obtained from (172|) by the substitution q^ 4-¥ l^. 

The second type of possible contributions to A M (g, I) can be obtained from (17UI) with 
the demand that the Lorentz index is carried by a matrix itself, 

7m> 757m, <V9"> v^l"- (73) 

The third type of terms from which the vertex A M (g, I) can be constructed contains 
the tensor e^ ual) 

e^XV, e^a^V, e^a^qpqT, e^a^l^P, q°P\, e„ 7 V^75- 

(74) 

Taking all terms (17TI) . (!72|) . (I731) and (!74|) together and using some 7 M algebra (for 
details see [HO]), it is possible to arrive to the most general expression for the vertex 

A M (g,0, 

A M (g, = /i(g 2 )g M + / 2 (g 2 )g M 7 5 + / 3 (<7 2 )7m + Z^h^s + Mq 2 )^ + / 6 (<? 2 ) W7^V, 

(75) 

where the only dependence on q 2 remains (because p 2 = p' 2 = m 2 where m is the fermion 
mass and I 2 = 4m 2 — q 2 ). 

From the natural requirement of current conservation (electromagnetic gauge invari- 
ance) d^j^ = it follows, that 

/i(?V + / 2 (g 2 )g 2 75 + 2m/ 4 (g 2 ) 7 5 = 0, (76) 

from which one gets 

/i(<7 2 ) = 0, f 2 (q 2 )q 2 + 2mf 4 (q 2 ) = 0. (77) 

Therefore, in the most general case consistent with Lorentz and electromagnetic gauge 
invariance, the vertex function is defined in terms of four form factors , 

Ml) = /o(? 2 )7m + f M (q 2 )i^q v + fjstf^fys + f A (q 2 )(q 2 l, - q„A)l5, (78) 

where fq(q 2 ), /m(<? 2 ), fE^q 2 ) and /a(<? 2 ) are charge, dipole magnetic and electric, and 
anapole neutrino form factors. 

Note that the form factors are Lorentz invariant and they depend only on q 2 , which 
is the only independent dynamical quantity which is Lorentz invariant. 

The hermiticity of the electromagnetic current and the assumption of its invariance 
under discreet symmetries transformations put certain constraints on neutrino form fac- 
tors, which are in general different for the Dirac and Majorana cases. In the case of Dirac 
neutrinos, the assumption of CP invariance combined with the hermiticity of the elec- 
tromagnetic current J^ M implies that the electric dipole form factor vanishes. At zero 
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Figure 14: Contribution of the neutrino vertex function to neutrino elastic scattering on 
a charged lepton. 



momentum transfer only /q(0) and /m(0), which are called electric charge and magnetic 
moments, contribute to the Hamiltonian, H int ~ J^A^, which describes the neutrino 
interaction with external electromagnetic field A M . It is also possible to show [5THI5T] that 
hermiticity by itself implies that /q, Jm, and Ja are real, 

Imf Q = Imfu = Imf A = 0. (79) 

In the case of Majorana neutrinos, regardless of whether CP-invariance is violated or 
not, the charge, dipole magnetic and electric form factors vanish [SUET], 

fq = fu = f E = 0. (80) 

This means that in the case of Majorana neutrinos only the anapole moment can be non- 
vanishing among the electromagnetic moments (see also [EE])- Note that it is possible 
to prove [57] that the existence of a non vanishing magnetic moment for a Majorana 
neutrino would bring a clear indication of CPT nonconservation. 

In general the matrix element of the electromagnetic current (1691) can be considered 
between different neutrino initial ipi(p) and final ipj{p') states of different masses, p 2 = 

2/2 2 
m i 1 V = m j- 

< ^(P')\J» M MP) >= ujWApMuib), (81) 
and the correspondent vertex function is defined in the most general form 

M<?) = (fQ(q 2 )ij + fA{q 2 )ijl^j {q 2 lt, ~q^A) + fhM 2 )^^ + fE(q 2 )ijO "75- (82) 

The form factors are matrices in the space of neutrino mass eigenstates [51]. General 
properties of the form factors in the diagonal case (i = j) have been already discussed. 
In the off-diagonal case (i 7^ j) the hermiticity by itself does not imply restrictions on 
the form factors of Dirac neutrinos. It is possible to show [57J that if the assumption 
of CP invariance is added, the form factors /q(<7 2 ), /m(<? 2 ), fE^q 2 ) and /a(9 2 ) should 
be relatively real to each other (no relative phases exist). For the Majorana neutrino, if 
CP invariance holds, there could be either a transition magnetic or a transition electric 
moment but not both. The anapole form factor of a Majorana neutrino can be nonzero. 
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3.1 Neutrino form factors in gauge models 

From the demand that the form factors at zero momentum transfer, q 2 = 0, are 
elements of the scattering matrix, it follows that in any consistent theoretical model 
the form factors in the matrix element (I6H|) should be gauge independent and finite. 
Then, the form factors values at q 2 = determine the static electromagnetic properties 
of the neutrino that can be probed or measured in the direct interaction with external 
electromagnetic fields. This is the case for charge, dipole magnetic and electric neutrino 
form factors in the minimally extended Standard Model . The neutrino anapole form 
factor is an exceptional case (see, for instance, [H2HS1]) arid will be discussed later in 
Section 13.41 

In non-Abelian gauge theories, the form factors in the matrix element fj69|) at nonzero 
momentum transfer, q 2 ^ 0, can be not invariant under the gauge transformation. This 
happens because in general the off-shell photon's propagator is gauge dependent. There- 
fore, the one-photon approximation is not enough to get physical quantities. In this case 
the form factors in the matrix element f[6T?j) cannot be directly measured in an experi- 
ment with an external electromagnetic field, however they can contribute to high order 
diagrams describing some processes that are accessible for experimental observation (for 
a discussion on this item see, for instance, [65]). As an example, a diagram for a neu- 
trino elastic scattering on a charged lepton is shown in Fig. where the hatched plaque 
represents the neutrino electromagnetic vertex function that includes contributions from 
the form factors. 

It should be noted that there is an important difference between the electromagnetic 
vertex function of massive and massless neutrinos [66] • For the case of a massless neutrino, 
the matrix element of the electromagnetic current (I6H|) can be expressed in terms of only 
one Dirac form factor fo{.q 2 ) (see, for example, also [SH]), 

u{p')^{q)u{p) = f D (q 2 Mp') lfl (l + l5 )u(p). (83) 

It follows that the electric charge and anapole form factors for a massless neutrino are 
related to the Dirac form factor /d(? 2 ) arid hence to each other 

f Q (q 2 ) = f D (q 2 ), f A (q 2 ) = f D (q 2 )/q 2 . (84) 

In the case of a massive neutrino, there is no such simple relation between electric 
charge and anapole form factors since the q^/rfb term in the anapole part of the vertex 
function (175|1 cannot be neglected. 

Consider [66] the full set of one-loop Feynman diagrams contributing to the Dirac 
massive neutrino electromagnetic vertex function in the framework of the Standard Model 
supplied with the S*?7(2)-singlet right-handed neutrino in the general gauge. The 
vertex function A At (g), in the one-loop approach, contains contributions given by two 
types of diagrams: the proper vertices (Fig. [To"]) and the 7 — Z self-energy diagrams 
(Fig. m. 

The direct calculation [66] of the massive neutrino electromagnetic vertex function, 
taking into account all of the diagrams (Fig. [15] and Fig. [T7|) . reveals that each of the 
Feynman diagrams gives nonzero contribution to the term proportional to 7^75. These 
contributions are not vanishing even at q 2 = 0. Therefore in addition to the usual four 



23 




X ' \ X 




(e) 



(f) 



Figure 15: (a) (f) Contributions to the neutrino vertex function from proper vertices (x is 



the unphysical would-be charged scalar boson; the correspondent Feynman rules necessary 
for the massive neutrino electromagnetic vertex calculations can be found in 
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Figure 16: Contributions to the neutrino vertex function of 7 — Z self-energy diagrams. 
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terms in (I78|) an extra term proportional to 7^75 appears and the corresponding additional 
form factor f§(q 2 ) can be introduced. This problem is related to the decomposition of the 
massive neutrino electromagnetic vertex function. The calculation of the contributions of 
the proper vertex diagrams (Fig.[T5|) and 7 — Z self-energy diagrams (Fig. fTB]) for arbitrary 

gauge fixing parameter a = 4 and arbitrary mass parameter a = 773- shows that at least 

s W 

in the zeroth and first orders of the expansion over the small neutrino mass parameter 
b = {j^) the corresponding "charge" (p = f$(q 2 = 0) is zero. The cancellation of 
contributions from the proper vertex and self-energy diagrams to the form factor f<o{q 2 ) 
at q 2 ^ 0, 

h(q 2 ) = ft Z) (Q 2 ) + ft° P - VCTt -\q 2 ) = 0, (85) 

was also shown [BB] for arbitrary mass parameters a and b in the 't Hooft-Feynman gauge 
a — 1. 

For a massive Dirac neutrino, by performing the direct calculations [66] of the complete 
set of one-loop diagrams it is established that the neutrino vertex function consists of only 
three electromagnetic form factors (in the case of a model with CP conservation). Closed 
integral expressions are found for electric, magnetic, and anapole form factors of a massive 
neutrino. On this basis, the electric charge (the value of the electric form factor at zero 
momentum transfer), magnetic moment, and anapole moment of a massive neutrino have 
been derived. It has been shown by means of direct calculations for the case of a massive 
neutrino that the electric charge is independent of the gauge parameters and is equal to 
zero, the magnetic moment is finite and does not depend on the choice of gauge. 



3.2 Neutrino electric charge 

It is usually believed [67] that the neutrino electric charge is zero. This is often thought 
to be attributed to gauge invariance and anomaly cancellation constraints imposed in the 
Standard Model . In the Standard Model of SU (2)l x U (l)y electroweak interactions it is 
possible to get [65|l6U| a general proof that neutrinos are electrically neutral. The electric 
charges of particles in this model are related to the SU{2)l and U{1)y eigenvalues by 
(see Table [T]) 

Q = h + \ (86) 

In the Standard Model without right-handed neutrinos z/r the triangle anomalies cancel- 
lation constraints (the requirement of renormalizability) lead to certain relations among 
particles hypercharges Y, that are enough to fix all Y, so that hypercharges, and conse- 
quently electric charges, are quantized [68J. In this case, neutrinos are electrically neutral. 

The direct calculation of the neutrino charge in the Standard Model under the assump- 
tion of a vanishing neutrino mass in different gauges and with use of different methods 
is presented in [6B* |[701 - [72] . For the flavor massive Dirac neutrino the one-loop contri- 
butions to the charge, in the context of the minimal extension of the Standard Model 
within the general gauge, were considered in [66] . By these direct calculations within 
the mentioned above theoretical frameworks it is proven that at least at one-loop level 
approximation neutrino electric charge is gauge independent and vanish. 

However, if the neutrino has a mass, the statement that a neutrino electric charge is 
zero is not so evident as it meets the eye. It is not entirely assured that the electric charge 
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should be quantized (see [73] and references therein). We recall here that the problem 
of charge quantization has been always a mystery within quantum electrodynamics [71] . 
The absence of an algebraic quantization of the charge eigenvalues in electrodynamics led 
to the proposal [75] of a possible topological explanation leading to magnetic monopoles. 

The strict requirements for charge quantization may also disappear in extensions of 
the standard SU(2)l x U(1)y electroweak interaction models if right-handed neutrinos 
ur with y ^ are included. In this case the uniqueness of particles hypercharges Y is 
lost (hypercharges are no more fixed) and in the absence of hypercharge quantization the 
electric charge gets "dequantized" [68]. As a result, neutrinos may become electrically 
millicharged particles. 

In general, the situation with charge quantization is different for Dirac and Majorana 
neutrinos. As it was shown in [69] . charge dequantization for Dirac neutrinos occurs in 
the extended Standard Model with right-handed neutrinos vr and also in a wide class 
of models that contain an explicit U(l) symmetry. On the contrary, if the neutrino is a 
Majorana particle, the arbitrariness of hypercharges in this kind of models is lost, leading 
to electric charge quantization and hence to neutrino neutrality [69J. 

Finally, while there are other Standard Model extensions (superstrings, GUT's etc) 
that provide enforcing of charge quantization, there are also models (for instance, with 
a "mirror sector" [76]) that predict the existence of new particles of arbitrary mass and 
small (unquantized) electric charge, in which neutrino can be a millicharged particle. 

The most severe experimental constraints on the electric charge of the neutrino 



are obtained assuming electric charge conservation in neutron beta decay n — > p-\-e~ + v e , 
from the neutrality of matter (from the measurements of the total charge q p + q e ) [77] 
and from the neutrality of the neutron itself [78]. Constraints from direct accelerator 
searches, charged leptons anomalous magnetic moments, stellar astrophysics and primor- 
dial nucleosynthesis are in general less stringent [7^179"]: 



A detailed discussion of different constraints on the neutrino electric charge can be found 
in [73]. 

3.3 Neutrino charge radius 

Even if the electric charge of a neutrino is vanishing, the electric form factor fq{q 2 ) can 
still contain nontrivial information about neutrino static properties. A neutral particle 
can be characterized by a superposition of two charge distributions of opposite signs 
so that the particle's form factor fq{q 2 ) can be non zero for q 2 ^ 0. The application 
of this notion to neutrinos has a long-standing history and is puzzling. In the case of 
an electrically neutral neutrino, one usually introduces the mean charge radius, which is 
determined by the second term in the expansion of the neutrino charge form factor fq{q 2 ) 
in series of powers of q 2 , 



q u < xlO-^e, 




q v < xlO 





(89) 
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The definition of the neutrino charge radius follows an analogy with the elastic electron 
scattering off a static spherically symmetric charged distribution of density p(r) (r = |x|), 
for which the differential cross section is determined [80H82] by the point particle cross 
section ^ 



«uJ /( ' 2)|2 - (90) 

where the correspondent form factor f(q 2 ) in the so-called Breit frame, in which go = 0, 
can be expressed as 

f(q 2 ) = J p(r)e iqx c/ 3 x = 4vr J drr 2 p(rf-^^, (91) 
here q = |q|. Thus, one has 

f - / r(r f™%^ r) ^. (92) 
In the case of small g, we have lim g 2^ qr cos ^ r 3 /"^ m ^^ = — y and 

f(q 2 ) = 1 - |q| 2 ^ + - • (93) 



Therefore, the neutrino charge radius (in fact, it is the charge radius squared) is usually 
defined by 

(rl) = -6^1,-0- (M) 

Since the neutrino charge density is not a positively defined quantity, (r%) can be negative. 

Just in one of the first studies [65], it was claimed that in the Standard Model and 
in the unitary gauge the neutrino charge radius is ultraviolet-divergent and so it is not 
a physical quantity. A recent direct one-loop calculation [66J of proper vertices (Fig. [T5|) 
and 7 — Z self-energy (Fig. ITB]) contributions to the neutrino charge radius performed in a 
general gauge for a massive Dirac neutrino gave also a divergent result. However, it was 
shown [83], using the unitary gauge, that by including in addition to the usual terms also 
contributions from diagrams of the neutrino-lepton neutral current scattering (Z boson 
diagrams), it is possible to obtain for the neutrino charge radius a gauge dependent but 
finite quantity. Later on it was also shown [39] that in order to define the neutrino 
charge radius as a physical quantity one has also to consider box diagrams (see Fig. [TSi) . 
which contribute to the scattering process v\ + I' — > v\ + and that in combination 
with contributions from the proper diagrams it is possible to obtain a finite and gauge- 
independent value for the neutrino charge radius. In this way, the neutrino electroweak 
radius was introduced [72] and an additional set of diagrams that give contribution to 
its value was discussed in [53]. Finally, in a series of recent papers [HI] the neutrino 
electroweak radius as a physical observable has been introduced. In the correspondent 
calculations, performed in the one-loop approximation including additional terms from 
the 7 — Z boson mixing and the box diagrams involving W and Z bosons, the following 
gauge-invariant result for the neutrino charge radius have been obtained: 
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(r 2 ) = 4 x l(T 33 cm 2 , (96) 



Figure 18: Contribution of the W box diagram to the scattering process vi + I' — > vi + V . 

(where mw and m ; are the W boson and lepton masses, % — e, [M, t). This result, however, 
revived the discussion [SoTISB] on the definition of the neutrino charge radius. Numerically, 
for the electron neutrino electroweak radius it yields 

"I ) = 4 x 10" 33 cm 2 

which is very close to the numerical estimation obtained much earlier in 

Note that the neutrino charge radius can be considered as an effective scale of the 
particle's "size", which should influence physical processes such as, for instance, neutrino 
scattering off electron (the differential cross section is given in Eq. (11211) below). To 
incorporate the neutrino charge radius contribution in the cross section, the following 
substitution [ST] can be used: 

1 2 
g v - + 2 sin 2 6 W + -m 2 w (rl e ) sin 2 9 W . (97) 

It is interesting to compare the theoretical results for the neutrino charge radius with 
some available experimental bounds [SB]: from primordial nucleosynthesis, 



(r 2 J < 7 x 10" 33 cm 2 , (98) 

from SN 1987A, 

(r 2 J < 2 x 10" 33 cm 2 , (99) 
from neutrino neutral-current reactions, 

- 2.74 x 10~ 32 cm 2 < (r 2 u J < 4.88 x 10" 33 cm 2 , (100) 
from solar experiments (Kamiokande II), 

(r 2 J < 2 x 10~ 32 cm 2 . (101) 

Recently, a new constraint have been obtained [BH] from a new evaluation of the weak 
mixing angle sin 2 9w by a combined fit of all electron neutrino elastic scattering data, 

- 1.3 x 10" 32 cm 2 < (r 2 e ) < 3.32 x 10~ 32 cm 2 . (102) 
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Comparing the theoretical value in Eq. ( 196]) with the experimental limits in Eqs. (198)) - 
f ll02p . one can see that they differ at most by one order of magnitude. Therefore, one 
may expect that the experimental accuracy will soon reach the value needed to probe the 
neutrino effective charge radius. 

It is obvious that the effects of new physics beyond the Standard Model can also 
contribute to the neutrino charge radius. In this concern, a recent work [59] should 
be mentioned, where the anomalous WW'y vertex contribution to the neutrino effective 
charge radius has been studied and the value for the correspondent additional contribution 
of 

|«)|<10- 34 cm 2 (103) 

was obtained. Note that this is only one order of magnitude lower than the expected 
value of the charge radius in the Standard Model . 

A detailed discussion on possibility to constrain the v T and charge radii from 
astrophysical and cosmological observations and from the terrestrial experiments can be 
found in [9"U] . 

3.4 Neutrino anapole moment 

The anapole form factor is the most mysterious and ambiguous among the neutrino 
form factors. The notion of an anapole moment for a Dirac particle was introduced 
in [HI] for a T-invariant interaction which, however, is not invariant under P and C 
transformations. 

To understand the physical meaning of the anapole form factor, as well as the meaning 
of other form factors, it is instructive to couple the correspondent term of the current to 
an external electromagnetic field (given by a potential A^), to derive the corresponding 
Dirac equation of motion for a neutrino field ip of mass m, and finally to obtain the 
interaction energy with a static electromagnetic field in the nonrelativistic limit. From 
this perspective, it is straightforward to understand that the charge form factor fq{q 2 ) 
at q 2 = is the electric charge, fq(q 2 ) = Q [821 [92]. Similarly, \i = /a/(0) and e = 
*/e(0) are t ne dipole magnetic and electric moments, respectively. In the nonrelativistic 
approximation, from the anapole term of the neutrino current (see Eqs. ([HUD and (1781) ). 
it is possible to obtain [51] the interaction energy 

H int K f A (0) (*■ curl B-E), (104) 

which corresponds to a T-invariant toroidal (anapole) interaction of the neutrino that does 
not conserve the P and C parities. This interaction defines the axial-vector interaction 
with an external electromagnetic field. The poloidal currents on a torus can be considered 
geometrical model for the anapole [93] . 
The direct calculation fSB] of the corresponding vertex contributions (the diagrams 
in Figs, [m and [T5]) to the massive Dirac neutrino anapole moment gives an infinite and 
gauge-dependent result. The same behavior of the charged leptons anapole moments has 
been demonstrated in [91] . Note that even in the case of massless neutrinos this is not a 
trivial task to obtain the anapole moment as a physical quantity (see section [373]) . Here 
we also recall that for the massless case the neutrino anapole moment is connected to the 
derivative of the electric charge form factor fq(q 2 ) with respect to q 2 at q 2 = 0, which is 
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the charge radius, by the relation 

av = f A (0) = Url). (105) 
o 

This relation is obtained within the Standard Model and in general it is model dependent. 
As it has been shown in [59J, the same relation between a massless Dirac neutrino anapole 
moment and charge radius in the context of an effective Yang-Mills theory which includes 
a general S'L r i / (2)-invariant Lorentz tensor structure of nonrenormalizable type for the 
WW'y vertex is also fulfilled. This relation is obtained within the Standard Model and in 
general it is model dependent. As it has been shown in [59], the same relation between a 
massless Dirac neutrino anapole moment and charge radius in the context of an effective 
Yang-Mills theory which includes a general S , 6 r i(2)-invariant Lorentz tensor structure of 
nonrenormalizable type for the WW'y vertex is also fulfilled. 

As it was discussed in jM], since the anapole form factor does not correspond to a 
multipole distribution, the anapole moment has a quite intricate classical analog. A more 
convenient and transparent characteristic, the toroidal dipole moment, was proposed in- 
stead for the description of T-invariant interactions. In this case, the electromagnetic 
vertex of a neutrino can be rewritten in an alternative multipole (toroidal) parameteri- 
zation. In some sense this parameterization has a more transparent and clear physical 
interpretation, because it provides a one-to-one correspondence between the multipole 
moments and the corresponding form factors. In one-loop calculations [M] of the toroidal 
(and anapole) moment of a massive and massless Majorana neutrino (the diagrams in 
Figs. [T5] and [THl contribute) it was shown that its value does not depend significantly on 

the neutrino mass (through the parameters ^r) and is of the order of 

m w 

Tu = fA(q 2 ) g i=o oc 10" 33 - 1CT 34 cm 2 , (106) 

depending on the values of the quark masses that propagate in the loop diagrams of 
Fig. El 

Note that the anapole form factors can contribute to the neutrino vertex function in 
both the diagonal and and off-diagonal cases. The anapole and the toroidal parameteri- 
zations coincide in the case when the current is diagonal on the neutrino initial and final 
masses. 

To conclude this section, it should be mentioned that the anapole interactions of a 
Majorana as well as a Dirac neutrino are expected to contribute to the total cross section 
of neutrino elastic scattering off electrons, quarks and nuclei. Due to the fact that the 
anapole interaction conserves helicity, its contribution to the cross section is similar to 
that of the neutrino charge radius. In principle, these contributions can be probed in 
low-energy scattering experiments in the future. 

3.5 Neutrino magnetic and electric dipole moments 

The neutrino dipole magnetic and electric form factors (and the corresponding mag- 
netic and electric dipole moments) are theoretically the most well studied and understood 
among the form factors. They also attract a reasonable attention from experimentalists, 
although the neutrino magnetic moment predicted in the Standard Model is proportional 
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to the neutrino mass and therefore is many orders of magnitude smaller than the present 
experimental limits obtained in terrestrial experiments. 

As it has been mentioned before, the first calculations of the neutrino dipole moments 
within a minimal extension of the Weinberg-Salam model (with nonzero neutrino mass 
and with a right-handed neutrino vr) were performed [1HHSD1IS2] by evaluating the ra- 
diative diagrams (a) and (d) shown in Fig. [15j The explicit evaluation of the one-loop 
contributions to the neutrino dipole moments in the leading approximation over the small 
parameters bi = jjfi- (here m 8 are the neutrino masses, i = 1,2,3), that in addition ex- 
actly accounts for the dependence on the small parameters a t = -S- (I = e, /i, r), yields, 
for Dirac neutrinos 



Mr 
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{ l±r ^) E ( 107 ) 



where 
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ai of 



f( ai ) = - 1 + 2 - 2 ' In a t . (108) 

4L 1-ai [l-aiY U ~ a U J 

All the charged lepton parameters ai are small. In the limit a; < 1 one has 

/(a,) « ^(l- leu). (109) 



2 V 2 

From Eqs. (11071) and (I109p . the diagonal magnetic moment of Dirac neutrinos are given 

by HBH5D] 

^^0^1^^). (110) 

1= e,[i,T 

Several important features of this result should be mentioned. The magnetic moment 
of a Dirac neutrino is proportional to the neutrino mass and for a massless Dirac neutrino 
in the Standard Model (in the absence of right-handed charged currents) the magnetic 
moment is zero. The magnetic moment of a massive Dirac neutrino, at the leading order 
in ai, is independent of the neutrino mixing matrix and also independent of the values of 
the charged lepton masses. The numerical value of the Dirac neutrino magnetic moment, 
as it follows from Eq. (jllOp . is 



^«3.2xl0- ia ^J/i B . (Ill) 

taking into account the existing constraints on neutrino masses, this value is several orders 
of magnitude smaller than the present experimental limits (see Section 13.61 for a further 
discussion on the experimental constraints on magnetic moments). 

From Eq. (II 07]) it can be clearly seen that in the Standard Model the static (diagonal) 
electric dipole moment of a Dirac neutrino vanishes, = 0. Dirac neutrinos may have 
nonzero diagonal electric moments in theories where CP invariance is violated. For a 
Majorana neutrino both the diagonal magnetic and electric moments are zero, ^ = 

e M = 0. 

'-it 

Let us discuss the neutrino transition moments, which are given by (11071) for i ^ 
j. If we again use the first two terms in the expansion (I109[) of the function f{a{) 
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and we insert the leading term in Eq. (I107j) . we get a vanishing result. This happens 
because the neutrino mixing matrix Uu is unitary and its rows and columns are orthogonal 
vectors. Therefore, the nonvanishing contribution comes only from the second term in 

the expansion of f(ai), which contains the additional small factor a/ = For the Dirac 
neutrino magnetic and electric transition moments, it is possible to obtain, rearranging 
the terms in Eq. fU07p . 





32 v / 2tt 2 V rn,iJ ^ \m w 

1= e, ^t, r 

Thus, they are reasonably suppressed with respect to the Dirac neutrino magnetic mo- 
ment (IllOp in the diagonal case (i = j). For convenience, numerically the Dirac transition 
moments can be expressed as follows (see, for instance, [73] ) 

4X10 - M ^) E ©W cm) 

1= e, r 

2 

The above-mentioned suppression by a factor of at least ai = -S- is due to the well- 
known Glashow-Iliopoulos-Maiani cancellation (GIM mechanism) [22] • Note that in the 
diagonal case (i = j) the leading term in the expression for the Dirac neutrino magnetic 
moment is not zero, because the sum in Eq. (11071) is equal to unity. 

Also Majorana neutrinos can have nonvanishing transition magnetic and electric mo- 
ments. In this case, additional Feynman diagrams should be considered, which also 
contribute to the dipole moments (for a detailed discussion see, for instance, [111153]). It 
is possible to show that, depending on the relative CP phase of the two neutrinos v\ and 
Uj, one of the two options is realized: fifj = 2/i^ and efj = 0, or nf- = and ef- = 2e®. 

In recent studies, the value of a massive Dirac neutrino diagonal magnetic moment 
was obtained in a one-loop approximation in the Standard Model , accounting for the 

2 

dependence on the neutrino mass parameter bi = 4^- [70] and accounting for the exact 

^ 2 

dependence on both mass parameters bi and ai = -jX- ES- The calculations of the 
neutrino magnetic moment which take into account exactly the dependence on the masses 
of all particles can be useful in the case of a heavy neutrino with a mass compared or even 
exceeding the values of other known particle masses. Note that the LEP data require 
that the number of light neutrinos coupled to the Z boson is exactly three. Therefore, any 
additional active neutrino must be heavier than ^4f-. In general, such a possibility is not 
excluded. That is the reason to consider the neutrino magnetic moment for various ranges 
of particles masses. The value of the neutrino magnetic moment for a light neutrino with 
mass m u « m f « M w that was obtained in [St^lTO]. 

&G 3 

H u = -=m u — --(2 — 7ai + 6af — 2af In a\ — af), (H4) 



4tt 2 V2 4(1 -at] 



reproduces the main term in Eq. (1 11 01) . i.e. the result derived in [4"5H5T)] . The authors of 
Ref. [66] obtained for an intermediate values of the neutrino mass, m u M\y, 

3eG F l\ 5 , 1 
H„ = ——^m„ { 1 + — b\ , (115) 



8tt 2 ^2 I 1 
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and for a heavy neutrino, rri£ <C M\y <c m v , 

_ 6__f (116) 

Note that in all the cases considered , the Dirac neutrino magnetic moment is proportional 
to the neutrino mass. This is an expected result, because the calculations have been 
performed within the Standard Model . 

In this concern, a question arises: "Is a neutrino magnetic moment always proportional 
to the neutrino mass?". The answer is "No". For example, much larger value for the 
Dirac neutrino magnetic moment can be obtained in SUl(2) x SUr{2) x U(l) left-right 
symmetric models (see, for instance, |HJ[96j[97] and the first paper in [71] ) with direct 
right-handed neutrino interactions. The intermediate gauge bosons mass states W\ and 
W2 have, respectively, predominant left-handed and right-handed coupling, since 

W 1 = W L coa£-W R wi\€ t (117) 

W 2 = W L sin£ + VP* cos f, (118) 
where £ here is a mixing angle and the fields Wl and Wr have pure V ± A interactions. 



The magnetic moment of a neutrino v% calculated in this model is 

m ,(! _ <L) sin2 £ + » (! (119 ) 

L v " tti/_ 7 t \ // III/. / J 



W2 t v ""W 2 ' 

where the term proportional to the charged lepton mass m/ is due to the left-right mixing. 
This term can exceed the second term in (11191) . which is proportional to the neutrino mass 



3.6 Experimental limits on neutrino magnetic moment 

The most sensitive and established method for the experimental investigation of the 
neutrino magnetic moment is provided by direct laboratory measurements of electron 
neutrino(antineutrino)-electron scattering at low energies in solar, accelerator and reactor 
experiments. A detailed description of different experiments can be found in [98|l99]. 

Extensive experimental studies of the neutrino magnetic moment, performed during 
many years, are stimulated by the hope to observe a value much larger than the prediction 
(II 111) of the minimally extended Standard Model (with nonzero neutrino masses). It 
would be a clear indication of new physics beyond the extended Standard Model . For 
example, in calculations [lOOj of the magnetic moment contribution to v e -e scattering in 
a class of extra-dimension models it was shown that the contribution to the cross section 
can be comparable with the corresponding one for the case in which there are no extra 
dimensions and the neutrino magnetic moment is of order \l v ~ lCT 10 /^. Future higher 
precision reactor experiments can therefore be used to provide new constraints on a class 
of large extra-dimension theories. 

The cross section for electron neutrino (antineutrino) scattering on electrons can be 
written |lUlj (see also [9^ |[9^l[TU2] ) as a sum of the Standard Model contribution and the 
neutrino magnetic moment contribution: 



da (da\ /da\ 



dt V dt 7 sm \dt 7 



(120) 
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Figure 19: Standard Model weak (W) and magnetic moment electromagnetic (EM) con- 
tributions to the cross section for several values of the neutrino magnetic moment [99J . 



The Standard Model contribution in Eq. ( I120p is 



dt J SM 



G 2 F m e 
2tt 



(gv + g A ) 2 + {gv - qa) 2 (i - -j^) 2 + (g 



2 m e T 



:i2i: 



where E v is the initial neutrino energy and T is the electron recoil energy, which is 
measured in experiments. The coupling constants gy and are 



gv 



2 sin 2 9 W + |, for v e , 



2 sin 6w 



2' 

i 

2 ' 



for i/ w ,i/ t 



9 a 



\, for i/ e , 
2, for v^v T . 



(122) 



In the case of antineutrinos, the substitution g^ —> —gA should be made. As it has been 
already mentioned, the neutrino charge radius can also contribute to the cross-section 
with the corresponding change of gy in Eq. 097|) . 

The neutrino magnetic moment contribution to the cross section is 



da 
~dt 



rat 



1-T/E U \ /^V 



(123) 



Note that the magnetic moment contribution to the cross section changes the helicity 
of the neutrino, contrary to the Standard Model contribution and also the possible con- 
tribution from the neutrino charge radius. Therefore, for relativistic neutrino energies 

is a negligible effect in the total cross see- 



the interference between ( % ) and ( % 

dt /SM \ dt s M 



tion f ll20p . The two terms an d [wj exm bit a quite a different dependence 

on the experimentally observable electron recoii energy T. The dependence of these two 
terms on T is shown [99] in Fig. [19] for six fixed values of the neutrino magnetic moment, 
= N x lCr n /iB, N = 1, 2, 3, 4, 5, 6. The cross sections are averaged over the typical 
antineutrino reactor spectrum (see also |lUlj ). It is easy to see that the lower the mea- 
sured recoil energy is, the smaller neutrino magnetic moment values are probed in the 
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experiment. From Eqs. ffT2T|) and (fT23j) . it follows that (^J exceeds (w) SM for 

T<^(^f- (124) 

The constraints on the neutrino magnetic moment in direct laboratory experiments 
have been obtained so far from the lack of any observable distortion of the recoil electron 
energy spectrum. Experiments of this type have started more than 30 years ago at the 
Savannah River Laboratory where the v-e scattering process was studied for the first 
time [103j . The upper limit on the magnetic moment fi u < (2-h4) x 10~ 10 //b was derived 
in Ref. |lUlj . The results of experiments at the Krasnoyarsk and Rovno reactors are, 
respectively, fi u < (2.4) x 10~ 10 fi B and /i u < (1.9) x KT 10 /^ [TU4] . The analysis of the 
recoil electron spectrum in the SuperKamiokande experiment gives \x v < (1.1) x 10 _10 /xb 
[105] . In reactor experiments carried recently, the following upper bounds have been 
obtained: ^ < 9.0 x 10- n /iB (MUNU [106]), ^ < 7.4 x 10~ u fi B (TEXONO [ID7j ). 
fi v < 5.8 x 10- u /ib (GEMMA pfh . The limit \i v < 5.4 x lO" 11 /^ has been recently 
obtained in the Borexino solar neutrino scattering experiment [109] . 

An upper limit on the neutrino magnetic moment \i v < 8.5 x 10~ n /iB has been found 
in an independent analysis of the first release of the Borexino experiment data performed 
in [HOj . It was also shown that with reasonable assumptions on the oscillation probability 
this limit translates into the upper limits on the magnetic moments of the muon and r 
neutrinos \i v < 1.5 x 10 _1( Vb and \i Vt < 1.9 x IO^Vb- The limit on u Ut is three order 
of magnitude stronger than that quoted by the Particle Data Group [lllj . 

An interesting new possibility for providing more stringent constraints on the neu- 
trino magnetic moment from v e -e scattering experiments was discussed in [112] on the 
basis of an observation |113j that the "dynamical zeros" appear in the Standard Model 
contribution to the scattering cross section. 

It should be mentioned that what is measured in experiments is an effective magnetic 
moment \x e ^ v whose value is a rather complicated function of the magnetic (transition) 
moments fiij. In addition, the dipole electric (transition) moments, if these quantities do 
not vanish, can also contribute to \f^ v . 

The magnetic moments fly, in the presence of mixing between different neutrino 
states, are associated with the neutrino mass eigenstates i^. The effective interaction 
Lagrangian which describes the coupling of neutrinos with the electromagnetic field is 
given by 

Lint = -4>iV a p{Hij + eijj^jF 01 ^ + h.c. (125) 

One can see that also the electric (transition) moments contribute to the coupling. 
In the laboratory neutrino scattering experiments on a neutrino magnetic moment the 
recoil energy of an electron, coupled to a neutrino flavor state, is the only measured 
quantity. That is why in order to extract from the experimental data information on 
the neutrino magnetic moment it is important to consider |114] the interplay between 
magnetic moment and mixing effects. Thus, the measured value fi exp depends on the 



The new stringent constraint on the level £t„ < 3.2 x 10 Vs has been also obtained within recently 
performed analysis |108j . 
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composition of the neutrino beam at the detector. The flavor composition of the initial 
beam of electron neutrinos is changing with the distance L from its source according to 



\v e (L)) = J2u ei e- iE * L Wi), (126) 



where E\ is the neutrino energy. Therefore the electromagnetic contribution to the scat- 
tering process amplitude is composed of terms like 



.1, - X! r "' " ; / 7<„- (127) 



In the magnetic scattering different mass eigenstates contribute incoherently, 

'da s 

' ii 



(I) „~£|£^- R v„ 
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Therefore, the effective value of the neutrino magnetic moment measured in scattering 
experiments is 

/4 p = plfa L, E v ) = E | E Uue'^^f. (129) 

3 i 

If a neutrino has nonvanishing electric (transition) moments, the substitution 

f-tij ->■ \fJ>ij - eij\ (130) 

should be made in Eq. (I129|) [73] . Therefore, in the case of Dirac neutrinos a destruc- 
tive interference between the magnetic and electric moments is possible. In the case of 
Majorana neutrinos, only magnetic or electric transition moments contribute to the cross 
section if CP is not violated. 

The general expression for fi^(ui,L,E v ) can be simplified |114] in several important 
cases. For instance, for Dirac neutrinos with only diagonal magnetic moments fiij = faS^ 
we have 

H 2 v (u e ,L,E v ) -> (/if) 2 = \ U eiW- (131) 

i 

Since in Eq. (I13ip there is no dependence on the distance L and the neutrino energy, 
the magnetic cross section is characterized by the initial neutrino flavor rather than 
by the composition of mass states in the detector. In this case, measurements of all 
"flavor" magnetic moments and mixing parameters, in principle, allow the extraction of 
the fundamental moments /ij. 

In the case of Majorana neutrinos, assuming that only two mass eigenstates are im- 
portant, we have 

ll? v (u e ,L,E v ) -> (/if) 2 = |/il 2 | 2 (|t/el| 2 + \Ue2\ 2 ) = |/il 2 | 2 , (132) 

which is independent on L, as well as on the neutrino energy and mixing. Note that the 
global fit |115] of the magnetic moment data from the reactor and solar neutrino exper- 
iments for the Majorana neutrinos produces limits on the neutrino transition moments 
A*23, A*si, A*ia < 1-8 x 10~ 10 /i B . 
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3.7 Model independent bounds on magnetic moments of Dirac 
and Majorana neutrinos 

As it was already mentioned before, there is a gap of many orders of magnitude 
between the present experimental limits ~ 10 _11 /i£ on neutrino magnetic moments (dis- 
cussed in the previous subsection) and the prediction f illip of the minimal extension of 
the Standard Model . At the same time, the experimental constraints have improved by 
only one order of magnitude during a period of about twenty years since the first limit 
on fi u was obtained [IQlj from neutrino scattering data. That is why experimental stud- 
ies of /i^ are in a reasonable extent stimulated by a hope that new physics beyond the 
minimally extended Standard Model might give much stronger contributions to fi u . One 
of the examples in which it is possible to avoid the neutrino magnetic moment being pro- 
portional to a (small) neutrino mass, that would in principle make a neutrino magnetic 
moment accessible for experimental observations, is realized in the left-right symmetric 
models considered before. 

Other interesting possibilities of obtaining neutrino magnetic moments lager than 
the prediction fillip of the minimal extension the Standard Model have been considered 
recently. In this concern, we note that it was proposed in |1UU] to probe a class of large 
extra dimensions models with future reactors searches for neutrino magnetic moments. 
The results obtained within the Minimal Supersymmetric Standard Model with i?-parity 
violating interactions |116] show that the Majorana transition magnetic moment might 
be significantly above the scale of f illip . 

Considering the generic problem with large neutrino magnetic moment, one can write 
down [55llll7j a naive relationship between the size of /i„ and the neutrino mass m u . 
Suppose that a large neutrino magnetic moment is generated by physics beyond a minimal 
extension of the Standard Model at an energy scale characterized by A. For a generic 
diagram corresponding to this contribution to fi u , one can again use the Feynman graph in 
Fig. EH the shaded circle in this case denotes effects of new physics beyond the Standard 
Model. The contribution of this diagram to the magnetic moment is 

fi v ~ (133) 

where e is the electric charge and G is a combination of coupling constants and loop 
factors. The same diagram of Fig. 13 but without the photon line gives a new physics 
contribution to the neutrino mass 

5m v ~ GA. (134) 
Combining the estimates f)133p and (11341) . one can get the relation 



between the one-loop contribution to the neutrino mass and the neutrino magnetic mo- 
ment. The A 2 dependence in Eq. f)135p was also discussed in |118U119j . 

It follows that, generally, in theoretical models that predict large values for the neu- 
trino magnetic moment, simultaneously large contributions to the neutrino mass arise. 
Therefore, a particular fine tuning is needed to get a large value for the neutrino mag- 
netic moment while keeping the neutrino mass within experimental bounds. One of the 
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possibilities |118] is based on the idea of suppressing the ratio m v j fi u with a symmetry: 
if a SU(2) U symmetry is an exact symmetry of the Lagrangian of a model, because of 
different symmetry properties of the mass and magnetic moment even a massless neu- 
trino can have a nonzero magnetic moment. If, as it happens in a realistic model, the 
SU(2) U symmetry is broken and if this breaking is small, the ratio m v / ' \x v is also small, 
giving a natural way to obtain a magnetic moment on the order of ~ 10 -11 /i£ without 
contradictions with neutrino mass experimental constraints. Several possibilities based 
on the general idea of |118] were considered in |120j . 

Another idea of neutrino mass suppression without suppression of the neutrino mag- 
netic moment was discussed in |119] within the Zee model |121j . which is based on the 
Standard Model gauge group SU(2) L x U(l)y and contains at least three Higgs doublets 
and a charged field which is a singlet of SU(2)l. For this kind of models there is a 
suppression of the neutrino mass diagram, while the magnetic moment diagram is not 
suppressed. 

It is possible to show with more general and rigorous considerations that the A 2 
dependence in Eq. (I135p arises from the quadratic divergence in the renormalization of the 
dimension-four neutrino mass operator. A general and model-independent upper bound 
on the Dirac neutrino magnetic moment, which can be generated by an effective theory 
beyond the Standard Model, has been derived [58] from the demand of absence of fine- 
tuning of effective operator coefficients and from the current experimental information 
on neutrino masses. A model with Dirac fermions, scalars and gauge bosons that is 
valid below the scale A and respects the Standard Model SU(2)l x t/(l)y symmetry 
was considered. Integrating out the physics above the scale A, the following effective 
Lagrangian that involves right-handed neutrinos vr, lepton isodoublets and the Higgs 
doublet can be obtained: 

^// = Ef^ n V) + ^-, (136) 

where \i is the renormalization scale, n > 4 denotes the operator dimension and j runs 
over independent operators of a given dimension. At n = 4 a neutrino mass arises from 
the operator O^f 1 = L&isr, where $ = ia^* ■ In addition, if the scale A is not extremely 
large with respect to the electroweak scale, then an important contribution to the neutrino 
mass can arise also from the higher dimension operators. At this point it is important 



to note that the combination of the n = 6 operators appearing in the Lagrangian (11 36 ft 
contains the magnetic moment operator ua^ u uF flL ' and also generates the contribution 
5m u to the neutrino mass (a detailed discussion of this item is given in Solving 
the renormalization group equation from the scale A to the electroweak scale, one finds 
that the contributions to the neutrino magnetic moment and to the neutrino mass are 
connected to each other by 

D 16a/2G F mjm u sin 4 9 W 
1 = 9^1/1 In (A/,) (137) 

where a is the fine structure constant, v is the vacuum expectation value of the Higgs 
doublet, 

/ = 1 - r - - tan 2 9 W - -(1 + r) tan 4 9 W , (138) 
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Figure 20: Schematic diagrams for neutrino-photon processes [731198]. 



r is a ratio of effective operator coefficients defined at the scale A and is of order unity 
without fine-tuning. If the neutrino magnetic moment is generated by new physics at a 
scale A ~ 1 TeV and the corresponding contribution to the neutrino mass is Sm u < 1 eV, 
then the bound [i v < 10 -14 /i£ can be obtained. This bound is several orders of magnitude 
stronger than the constraints from reactor and solar neutrino scattering experiments 
discussed before. 

The model-independent limit on a Majorana neutrino transition magnetic moment was 
also discussed in [58J. However, the limit in the Majorana case is much weaker than that in 
the Dirac case, because for a Majorana neutrino the magnetic moment contribution to the 
mass is Yukawa suppressed. The limit on ii^f is also weaker than the present experimental 
limits if /A^f is generated by new physics at the scale A ~ 1 TeV. An important conclusion 
of [58], based on model-independent considerations of the contributions to //„, is that if 
a neutrino magnetic moment of order \x v > 10~ 15 /iB were observed in an experiment, it 
would give a proof that neutrinos are Majorana rather than Dirac particles. 

4 Effects of neutrino electromagnetic properties 

If a neutrino has non-trivial electromagnetic properties, notably nonvanishing mag- 
netic and electric (transition) dipole moments or nonzero millicharge and charge radius, 
then a direct neutrino coupling to photons is possible and several processes important 
for applications exist. A set of typical and most important neutrino electromagnetic pro- 
cesses involving the direct neutrino couplings with photons is shown in Fig. [201 These 
processes are: 1) the neutrino radiative decay V\ — > ^2 + 7; neutrino Cherenkov radia- 
tion in an external environment (plasma and/or electromagnetic fields), the spin light 
of neutrino, SLu, in the presence of a medium; 2) the photon (plasmon) decay to a 
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neutrino- antineutrino pair in a plasma 7 — > vv\ 3) neutrino scattering off electrons (or 
nuclei), this process has been already considered in detail in Section [3T6| 4) neutrino spin 
(spin-flavor) precession in a magnetic field. 

4.1 Neutrino radiative decay and other v — > v + 7 processes 

If the masses of neutrinos are not degenerate, the radiative decay of a heavier neutrino 
v i into a lighter neutrino Vj (rrii > rrij ) with emission of a photon, 

Ui Uj + 7, (139) 

may proceed in vacuum [41 ^ |41? | I51 ^ I51? ] I5T) ]["1221I123] . A discussion of the possible role of 
the neutrino radiative decay in different astrophysical and cosmological setting has been 
started in [121]. 

The neutrino radiative decay process is described by the effective Lagrangian given in 
Eq. f)125p . The corresponding typical one-loop Feynman diagrams that contribute to the 
process in the standard SU(2)l x U(1)y model are similar to those shown in Fig. 15 if to 
consider the initial and final neutrinos as different mass states. For the case of a Dirac 
neutrino, the decay rate is found to be [^l3mi52ll531l551 fT22lfr23] 

aG% /mf — m?\3 1 2 

-f +7 = 12^4 (-^) K 2 + M E f^W^\ , (140) 



where /(a/) is given by Eq. (jl08p . Recalling the results for the Dirac neutrino magnetic 
and electric transition moments fiij and e^, given in Eq. (I107j) . one may rewrite Eq. (11 401) 
in the following form (see, for instance, [73]): 



kil 2 + l4l ( m \- m ) 



sir \ rrij 



(141) 



For degenerate neutrino masses (m 8 = rrij), the process is kinematically forbidden in 
vacuum. 

Note that there are models ( see for instance [125] ) in which the neutrino radiative 
decay rate (as well as the magnetic moment discussed above) of a non-standard Dirac 
neutrino are much larger than those predicted in minimally extended Standard Model. 

In the evaluation of the decay rate for Majorana neutrinos, two cases should be consid- 
ered that correspond to the two possible relative CP phases of vi and Vy If the Majorana 
neutrinos Vi and Vj have the same CP eigenvalues, the decay rate is 

r ^+T=6i^(^-^) K 2 -™?)| E f^ U i\- ( 142 ) 

^ 1= e, fi, t 

In this case the decay process is induced purely by the neutrino electric transition dipole 
moment, because ^ = 0. If the Majorana neutrinos have opposite CP eigenvalues, the 
decay rate is 



v-^ = ^( ! ^^) s (^ a +^)| E /(«otw| a - ( 143 ) 

1= e, fi, 



aG% /mf - mi 
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Figure 21: Astrophysical limits on neutrino transition moments [73] . 



In this case, the transition is purely of magnetic dipole type (e^ = 0). 



For numerical estimations it is convenient to express Eq. (I14ip in the following form: 



5.3 x 



2 /mi 



m„ 



m,- 



4 eV' 1 



(144) 



with the effective neutrino magnetic moment /i e // = w |yUjj| 2 + |e?-|. The expression (I144p 

is valid for both Dirac and Majorana neutrinos. In the case of Majorana neutrinos, only 
fiij or €ij contributes, depending on the relative CP phase of the neutrino states. Note 
that there is no destructive interference between and in Eq. (I144p . contrary to the 
neutrino scattering off electrons cross section given in Eqs. (j!23j) and (11301) . 

The neutrino radiative decay can be constrained by the absence of decay photons in 
reactor v e and solar v e fluxes. The limits on /x e // that are obtained from these consider- 
ations are much weaker than those obtained from neutrino scattering terrestrial exper- 
iments. Stronger constraints on fi e ff (though still weaker than those mentioned above) 
are obtained from the neutrino decay limit set by SN 1987A and also from the limit on 
the cosmic microwave background radiation distortions. These limits can be expressed 
as (see [73] and references therein) 



O.QxlO-H— V 
0.5xKT 5 (^y 
1.5xl(T 8 (^y 
1.0xl0- u (^) 9/4 



Reactor (z/ e ), 
Sun (u e ), 

SN 1987A (all flavors), 

Cosmic background (all flavors). 



(145) 



A detailed discussion (and corresponding references) on astrophysical constraints on the 
neutrino magnetic and electric transition moments, summarized in Fig. 21, can be found 
in [73]. 

For completeness, we would like to mention that other processes characterized by the 
same signature of Eq. (I139p have been considered previously (for a review of the literature 
see (I2SHI25] ): 

i) the photon radiation by a massless neutrino (z^ — > i/j +7, % — j) due to the vacuum 
polarization loop diagram in the presence of an external magnetic field |130j : 

ii) the photon radiation by a massive neutrino with nonvanishing magnetic moment 
in constant magnetic and electromagnetic wave fields [131] ; 
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iii) the Cherenkov radiation due to the nonvanishing neutrino magnetic moment in 
an homogeneous and infinitely extended medium, which is only possible if the speed of 
the neutrino is larger than the speed of light in the medium |132] ; 

iv) the transition radiation due to a nonvanishing neutrino magnetic moment which 
would be produced when the neutrino crosses the interface of two media with different 
refractive indices |133j ; 

v) the Cherenkov radiation of a massless neutrino due to its induced charge in a 
medium |134^l 

vi) the Cherenkov radiation of massive and massless neutrinos in a magnetized medium 
[125111571 ■ 

vii) the neutrino radiative decay (z^ — > Uj + 7, i ^ j) in external fields and media 
(see |138] and references therein). 

Recently, another mechanism of electromagnetic radiation by a massive neutrino in 
presence of matter (termed the spin light of neutrino, SLu), has been proposed |127j . 
The SLu is an electromagnetic radiation that can be emitted by a massive neutrino due 
to the neutrino magnetic or electric (transition) moments when the particle moves in the 
background matter. Within a quasi-classical treatment, the existence of the SLu was 
first studied |127[I139] on the basis of the developed Lorentz invariant approach to the 
neutrino spin evolution that implies the use of the generalized Bargmann-Michel-Telegdi 
equation |14Uj . 

Within the developed Lorentz invariant approach, it is also possible to find the solution 
of the neutrino spin evolution problem for a general case when the neutrino is subjected 
to general types of non- derivative interactions with external fields |141] (see also |142j ). 
These interactions are given by the Lagrangian 

-C = g s s{x)vv+g p ^{x)v 1 b v+g v V^{x)v 1 ^g a A^{x)v 1 

(146) 

where s,ir, = (V°,V),A fl = (A°,A),T^ = (a, b),n^ = (c, d) are the scalar, pseu- 
doscalar, vector, axial-vector, tensor and pseudotensor fields, respectively. For the corre- 
sponding spin evolution equation it has been found 

at I 1 + 7 1 7 J 

+2g t |[S x b] - + [S x [a x /3]]} (147) 

+%;j[Sxc]-^M-[Sx[dx/3]]}. 

This is a rather general equation for the neutrino spin evolution that can be also used for 
the description of neutrino spin oscillations in different environments, such as moving and 
polarized matter with external electromagnetic fields (see |128] and references therein). 
The SLu in gravitational fields has been studied (see the second paper of [139j ). for the 

3 Note that the neutrino electromagnetic properties are in general affected by the external environment. 
In particular, a neutrino can acquire an electric charge in magnetized matter |134j and the neutrino 
magnetic moment depends on the strength of external electromagnetic fields |135j . A recent study of the 
neutrino electromagnetic vertex in magnetized matter can be found in |136j . See also |128j for a review 
of neutrino interactions in external electromagnetic fields. 
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first time, on the basis of a neutrino spin evolution equation (11471) . Some general aspects 
of chiral dynamics for a neutrino non-minimally coupled with an external magnetic field 
have been discussed in |143j . 

It should be mentioned that the SLv in matter is a new mechanism of electromag- 
netic radiation that cannot be considered as the neutrino Cherenkov radiation in matter 
mentioned above, because it can exist even when the emitted photon refractive index is 
equal to unity. The SLv radiation is due to radiation of the neutrino by its own, rather 
than radiation of the background particles. As it was clear from the very beginning |127j . 
the SLv is a quantum phenomenon by its nature. The quantum theory of this radiation 
has been elaborated [144] (see also [145j ) within a development |129[I146[I147[I15U] of a 
quite powerful method that implies the use of the exact solutions of the modified Dirac 
equation for the neutrino wave function in matter. The corresponding Feynman diagram 
of the SLv processes is shown in Fig. [221 where the neutrino initial (ipi) and final ("0/) 




Figure 22: The spin light of neutrino (SLv) radiation diagram. 



states (shown by "broad lines" ) are exact solutions of the corresponding Dirac equations 
accounting exactly for the interaction with matter. The amplitude of the SLv process is 
given by 

/ ikx 
d 4 xjj f {x){te*) : ijjjjx), f = iu{ [S x x] + *7 5 £}, (148) 
V 2a; L 3 

where \x is the neutrino magnetic moment, = (u, k) and e* are the photon momentum 
and polarization vectors, x = k/cj is the unit vector pointing in the possible direction 
of the emitted photon propagation. From Eq. (1148)) it follows that the SLv rate and 
radiation power are proportional to /x^, so that these quantities are in general quite small. 
At the same time, for a wide range of matter densities these characteristics of radiation are 
increasing with the neutrino momentum. One may expect [127H1291H39pi44H150j that this 
radiation can be produced by high-energy neutrinos propagating in different astrophysical 
and cosmological environments. In the most interesting for possible astrophysical and 
cosmology applications case of ultra-high energy neutrinos, the average energy of the SLv 
photons is one third of the neutrino momentum, so that in principle the SLv spectrum 
spans up to the range peculiar of gamma-rays. 

It should be emphasized that the SLv mechanism of radiation (i.e. the transition 
between neutrino states with equal masses) can only become possible because of an 
external environment (plasma) influence on neutrino states. A possible impact of the 
background plasma on the SLv radiation through the plasma influence on propagation 
of SLv photons has been first considered in [144] . The plasma effects for the SLv were 
further studied in |151j where the role of the SLv plasmon mass was discussed. In the 



44 



case of ultra-high energy neutrino (i.e., in the only case when the time scale of the process 
can be much less than the age of the Universe) the SLv rate of |151] exactly reproduces 
the result obtained in [144] . For a more detailed discussion on the historical aspects of 
this issue see [T47HI41?] . 



4.2 Photon (plasmon) decay into neutrino-antineutrino pair 

The most interesting process, for the purpose of constraining neutrino electromagnetic 
properties, is the photon (plasmon) decay into a neutrino-antineutrino pair, 7* — > v + v. 
This plasmon process becomes kinematically allowed in media, because a photon with 
the dispersion relation u 2 + > roughly behaves as a particle with an effective mass. 
Note that the 7* — > v + v generated by the neutrino coupling to photons due to a 
magnetic moment fi u (and/or also due to a neutrino electric millicharge q u ) was first 
considered in [67] as a new energy-loss channels of the Sun. From the requirement that 
new energy-loss channels do not reasonably exceed the standard solar model luminosity, 
the constraints are found [73] to be \i v < 4 x 10~ 10 fi B and q u < 6 x 10~ 14 e. 

The tightest astrophysical bound on a neutrino magnetic moment is provided by 
observed properties of globular cluster stars. The plasmon decay 7* — > v + v inside 
the star liberates the energy w 7 in the form of neutrinos that freely escape the stellar 
environment. This nonstandard energy loss cools a red giant star so fast that it can delay 
helium ignition. The energy-loss rate per unit volume due to the plasmon decay process 
is 

9 



, tuyfk^^pd fey, (149) 

where is the photon Bose-Einstein distribution function, g — 2 is the number of 
polarization states and the decay rate is 

r 7 _.iiKzS! (150) 

24tt u; 7 

For a sufficiently large neutrino magnetic moment, the plasmon decay rate can be 
enhanced, inducing a significant delay of helium ignition. From the lack of observational 
evidence of this effect, the following limit has been found [152] : 

|i„<3x 1Q- 12 Hb- (151) 

This is the most stringent astrophysical constraint on a neutrino magnetic moment, ap- 
plicable to both Dirac and Majorana neutrinos. The same limit applies for the neutrino 
magnetic transition moments as well as for the electric (transition) moments. 



4.3 Neutrino spin-flavor precession in magnetic field 

If neutrinos have magnetic moments, the spin can precess in a transverse magnetic field 
[153H155] . Considering for simplicity only one undetermined flavor, a neutrino produced 
at x = is described at a distance x by the state 

W(x)) = <p L (x) \v l ) + <Pr(x) \vr) , (152) 
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where \vl) and \ur) are, respectively, neutrino states with negative and positive helicity, 
which are called left-handed and right-handed. The respective amplitudes ^>l{x) and 
<Pr(x) should not be confused with chiral fields. Their evolution equation in a transverse 
magnetic field B±(x) is given by 



dx 



M x ) 
M x ) 





fxB ± (x) 



fj,B ± (x] 




M x ) 
M x ) 



(153) 



where fi is the magnetic moment (we consider a Dirac neutrino). This differential equation 
can be solved through the transformation 



M x ) 
Vr{ x ) 



1 

71 



1 1 
-1 1 



(154) 



The new amplitudes <f-(x) and (p+(x) satisfy decoupled differential equations, whose 
solutions are 



f^yx) = exp 

If we consider an initial left-handed neutrino, we have 



fX 

±i / dx'fiB ± (x') 
Jo 



Mo)- 



Mo) 
Mo) 



Mo) 
Mo) 



Then, the probability of vl — > Vr transitions is given by 



\<Pr( x ) 



sm 



dx' \x B\ 



[x 



(155) 



(156) 



(157) 



Note that the transition probability is independent from the neutrino energy (contrary 
to the case of flavor oscillations) and the amplitude of the oscillating probability is unity. 
Hence, when the argument of the sine is equal to tt/2 there is complete vl — > vr conver- 
sion. 

The precession v eh — > u eR in the magnetic field of the Sun was considered in 1971 |153] 
as a possible solution of the solar neutrino problem. If neutrinos are Dirac particles, right- 
handed neutrinos are sterile and a v eL — > u eR conversion could explain the disappearance 
of active solar z/ eL 's. 



In 1986 it was realized |154[I155] that the matter effect during neutrino propagation 
inside of the sun suppresses v eL — > u eR transition by lifting the degeneracy of v eL and u eR . 
Indeed, taking into account matter effects, the evolution equation f)153p becomes 



M x ) 
dx Wr(x) 



V nB ± (x] 
fiB ± (x) 



M x ) 



(158) 



with the appropriate potential V which depends on the neutrino flavor, according to 
Eq. fj44l) . Again, we consider a Dirac neutrino, which can have a magnetic moment. In 
the case of a constant matter density, this differential equation can be solved analytically 
with the orthogonal transformation 



M x ) 
Vr( x ) 



cos£ sin£\ / <P-{x] 
-sin£ cosU \<f+(x] 



(159) 
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The angle £ is chosen in order to diagonalize the matrix operator in Eq. ( II 5 8 p : 



sin 2£ 



with the effective energy splitting in matter 



2/iB A 
AE- 



M 



AE 



M 



V 2 + (2/zSj 



The decoupled evolution of (p T (x) is given by 



<p=f(x) = exp 



MO) • 



(160) 



(161) 



(162) 



For an initial left-handed neutrino, 



Mo) 

MO) 



cos A 
sin£ ) ' 



(163) 



leading to the oscillatory transition probability 



P VL -+ VR {x) = \vr(x)\ = sin 2£sin ( -AE M x ) . 



(164) 



Since in matter AEm > 2fiB±, the matter effect suppresses the amplitude of vl ^r 
transitions. However, these transitions are still independent from the neutrino energy, 
which does not enter in the evolution equation 



Once it was known, in 1986 |154|I155] , that the matter potential has the effect of 
suppressing ul —> vr transitions because it breaks the degeneracy of left-handed and 
right-handed states, it did not take long to realize, in 1988 |156U157j . that the mat- 
ter potentials can cause resonant spin-flavor precession if different flavor neutrinos have 
transition magnetic moments (spin-flavor precession in vacuum was previously discussed 
in PI)- 

A detailed discussion on how to attack the solar neutrino problem using a neutrino 



magnetic moment can be found in |159U160j . 

Let us consider two neutrino flavors: v e and v^. A neutrino produced at x = is 
described at a distance x by the state 



\v{x)) = lf eL {x) \v eL ) +¥e R {x) We R ) + <P» L (x) \Vp L ) + <P m (x) |z^J 



(165) 



which is the generalization of Eqs. (T4"7"|) and (I152|) . 

Considering Dirac neutrinos, which can have diagonal magnetic moments /i ee and 
as well as transition magnetic moments fi efl and /i^, the evolution equation of the 



amplitudes, obtained by combining Eqs. 
terms, is 



and (11581) and neglecting irrelevant diagonal 



da; 



( <Pe L ( X )\ 
<fe R {x) 



( Ve L { x )\ 

<Phl( x ) 
<Pe R {x) 

\<Phr( x )J 



(166) 
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with the effective Hamiltonian matrix 



n 



f-^cos2d + V e ^sin2tf VeeB ± (x) ^B ± (x)\ 
A|i sin 2 ^ ^ cos 2d + V„ ^ e B ± {x) fi^B ± (x) 



fi ee B ± (x) 
\ ^B±{x) 



HfieB±(x) 



Am2 cos2^ 4#sin2tf 



4 

^ sin 2t9 ^fcos2dJ 



45 



(167) 



4£ 



The matter potential can generate resonances, when two diagonal elements of % become 
equal. There are two resonances: 



1 . There is a resonance in the u eL u IJiR channel for 

K = — — cos 2$. 
2£ 



(168) 



The density at which this resonance occurs is not the same as that of the MSW 
resonance, given by Eq. ( T56|) . because of the neutral-current contribution to V e = 
Vcc + Vnc- The location of this resonance depends on both N e and N n . 



2. There is a resonance in the v, 



I'L 



v eR channel for 



Am 2 
'~2E 



cos 2d. 



(169) 



If cos 2d > 0, this resonance is possible in normal matter, since the sign of = Vnc 
is negative, as one can see from Eq. (l4"5l) . 



In practice the effect of these resonances could be the disappearance of active v eL or v iiL 
into sterile right-handed states. 

Let us consider now the more interesting case of Majorana neutrinos, which presents 
two fundamental differences with respect to the Dirac case: 

(A) Majorana neutrinos can have only a transition magnetic moment fi efl = —fi^ e - 

(B) The right-handed states are not sterile, but interact as right-handed Dirac antineu- 
trinos. 



The evolution equation of the amplitudes is given by Eq. f )166p with the effective Hamil- 
tonian matrix 



n 



(- ^ cos 2d + V e ^£ sin 2d 
^sin2tf ^g-cos2d + V ll 







^ HetiB±(x) 



-pi efl B ± (x) 




Again, there are two resonances: 







-^B ± (x) 

' COS 2d - V e 



Am 2 
4E u,B - iu 



^sin2^ 



fi ef ,B ± (x) \ 




Am 2 
AE 



sin 2d 



^cos2d-V, 



(170) 



I. There is a resonance in the u eL ^ v flR channel for 



^cc + 2V 



NC 



Am 2 
~2E 



cos 2d . 



(171) 
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2. There is a resonance in the v fXL £5 z/ efl channel for 

V^cc + 2^ NC = -— cos 2^ . (172) 

The location of both resonances depend on both N e and N n . If cos 2$ > 0, only the first 
resonance can occur in normal matter, where N n ~ N e /6. A realization of the second 
resonance requires a large neutron number density, as that in a neutron star. 

The neutrino spin oscillations in a transverse magnetic field with a possible rotation 
of the field-strength vector in a plane orthogonal to the neutrino-propagation direction 
(such rotating fields may exist in the convective zone of the Sun) have been considered in 
|158[ll61j . The effect of the magnetic-field rotation may substantially shift the resonance 
point of neutrino oscillations. Neutrino spin oscillations in electromagnetic fields with 
other different configurations, including a longitudinal magnetic field and the field of an 
electromagnetic wave, were first examined in |162] and |14Ul [163j. 

It is possible to formulate a criterion |158] for finding out if the neutrino spin (spin- 
flavor) precession is significant for given neutrino and background medium properties. 
The probability of oscillatory transitions between two neutrino states v a L ^ v^r can be 
expressed in terms of the elements of the effective Hamiltonian matrices (I167P and (I170p 
as ^ 

p v aj ^v, = sin 2 e// sin 2 , (173) 



eff 



where 



sin 2 & e ff = = ; — — — -, and L e rt = — , (174) 



The transition probability can be of order unity if the following two conditions hold 
simultaneously: 1) the amplitude of the transition probability must be "far" from zero 
(at least sin 2 $ e ^ > 1/2), 2) the neutrino path length in a medium with a magnetic field 
should be longer than half the effective length of oscillations L e ff. In accordance with 
this criterion, it is possible to introduce the critical strength of a magnetic field B cr which 
determines the region of field values B± > B cr at which the probability amplitude is not 
small (sin 2 ^ e// > 1/2): 

1 



Bar = —^(Hw-H aa ) 2 , (175) 

where fi is fj, ee , /i eAt , or /x Me depending on the type of neutrino transition process in 
question. 

Consider, for instance, the case of v eh ^ u flR transitions between Majorana neutrinos. 
From Eqs. flT75|) and (TT70|) . it follows [T5B] that 



B r 



U^A-V2G F N eff 
2/2 V 2E n 



;i76) 



where A = cos 2$ and N e ff = N e — N n . For getting numerical estimates of B cr it is 
convenient to re-write Eq. (11761) in the following form: 



B cr « 43^ 



leV 2 /V E u J 1 cm" 3 



Gauss. (177) 
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An interesting feature of the evolution equation (I166p in the case of Majorana neutri- 
nos is that the interplay of spin precession and flavor oscillations can generate v eh — > v eR 
transitions |164j . Since v eR interacts as right-handed Dirac antineutrinos, it is often de- 
noted by u eR , or only z/ e , and called "electron antineutrino" . This state can be detected 
in through the inverse /3-decay reaction 

v e +p^n + e + , (178) 

having a threshold E t h = 1.8 MeV. 

The possibility of v eh — > v eR transitions generated by spin-flavor precession is particu- 
larly interesting for solar neutrinos, which experience matter effects in the interior of the 
Sun in the presence of the solar magnetic field. In 2002, the Super-Kamiokande Collab- 
oration established for the flux of solar z/ e 's an upper limit of 0.8%, at 90% C.L., of the 
Standard Solar Model neutrino flux in the range of energy from 8 to 20 MeV [165] . This 
limit was improved in 2003 by the KamLAND Collaboration to 2.8 x 10 -4 in the energy 
range 8.3 - 14.8 MeV [166] . The implications of this limit for the spin-flavor precession 
of solar neutrinos have been studied in several papers [167H171] , taking into account 
the dominant v e — > u^, v T transitions due to neutrino oscillations (see the brief review 
in Section [22]). Considering turbulent solar magnetic field models in which v e ^ — \ v e ^ 
transitions are strongly enhanced, the authors of Refs. [167U168] obtained the interesting 
limit 

H ea < few x 10^ 12 At B , (179) 

where \x ea is the transition magnetic moment between v e and v a = cos $23^ — sin #23 ^r- 
This limit has been, however, criticized in Ref. [170] . 

The spin-flavor mechanism was considered [172] in order to describe time variations of 
solar-neutrino fluxes in gallium experiments. The effect of a nonzero neutrino magnetic 
moment is also of interest in connection with the analysis of helioseismological observa- 
tions [173] . 

The idea that the neutrino magnetic moment may solve the supernova problem, i.e. 
that the neutrino spin-flip transitions in a magnetic field provide an efficient mechanism of 
energy transfer from a protoneutron star, was first discussed in [174] and then investigated 
[175] in some detail. The possibility of a loss of up to half of the active left-handed 
neutrinos because of their transition to sterile right-handed neutrinos in strong magnetic 
fields at the boundary of the neutron star (the so-called boundary effect) was considered 
in [T58] . 

In conclusion, we would like to point out [176] that there is a huge gap of many orders 
of magnitude between the present limits oc lO^ 11 ^ 14 )/!^ on a neutrino magnetic moment 
[i v and the prediction of a minimal extension of the Standard Model. Therefore, if any 
direct experimental confirmation of non-zero neutrino magnetic moment were obtained 
within a reasonable time in the future, it would open a window to new physics. 
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